Analysis of Yale Cold-Fusion Experiment 
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Introduction , , , . . * 

The following represents an analysis of some of the data from a cold-fus,on expenmen 
A YrIp in the fall of 1989. The purpose of the experiment was to see if proton-recoi 

coutd'p^de adequate information about neutron 

singles rates, or the time-structure and size of the neutron bursts which ave een 

in other experiments. 


The Experiment . 

The data consist of puls, heights and puW shape, (vridth.) for tao ggwW 
of neutron detectors, e«h consisting rf mre central Z 

ring counters of ^8% f < 

wheTCThat detector had fired. The trigger requirement was apparently a hit in either 

This analysis has been made difficult by the lack of any logic or 
for the experiment, which necessitated mud. detective work and frequent d.gress.ons. 


Neutron Singles .♦ . . , ... „„ 

With moderately -constricting twcedimensional cuts 

sus pulse shape, unvetoed neutron-ringlee rate, were the orf« of lt^O/to tm ea« 
debtor for both foreground and background runs. Putting t!ght time^-flight ^ts on 
center-ring coincidences provided a measure of neutron efficiency and yielded o 

cuts Th<f quality of pulse-shape discrimination between neutrons and gammas was very 
^ecTX at low pulse heights, and in both data and background runs, very few of 

the surviving events appeared to be gammas. . , •.«(* n t AUg 

Although the veto counters were advertised as only 65% efficient, only 10% 
single events were vetoed, indicating that nuddng the veto more effiaent wril probably 

not gain much. 
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• f ted nulse-shape/pulse-height spectra for background 

Examination of gated l p^e shap / > reol neu trons. r 

Shows that !%+%*£££& le. «** 
detectors yields the toliowing u imvNp1lt 

tvoe 

Table 1: 


background runs (T'S L ') 
Looking just at the center 


run 

42 

43 

44 

45 

46 

47 
53 

total 


type 

FG 

FG 

FG 

FG 

FG 

FG 

FG 

FG 


lOOxNeut/Good 

4.92 ±0.12 
5.02 ±0.12 
4.78 ±0.12 
4.83 ±0.12 
4.95 ±0.12 
4.60 ±0.12 
4.75 ±0.12 
4.86±0.04 


40 

41 
53 
total 


BG 

BG 


4.88 ±0.12 
5.08 ±0.15 
nn 4.53 ±0.16 

w BG 4.84 ±0.10 

Note the variation cd as much as 

measurements. ancles rates, and the variability of background, it is neces ary 

to iScidlt can reduce the bakgo md -^ ^ „f coincidence. ehoUd be 

Note that for real neutron, of P*?P“ erLta The only purpose coincident 

•— * e — : 3 

s value x , , 

< 0 Bad pulse 

^ Neutron 

o Gamma 

2 0K both latcnes ^ , 

The P1D value in the graph i. then P1D = ^toTlHLarared 

events are con^^d » D5 (#10) TOFia lO^hUe that for U3 is 

Wow « TDC proems 1" "era 3E ^ "te 

* Av^ndix A is similar to Table 1, but tor nn S . background neutrons to 

tzszrjz-**.- a——- 


CR 

no latches 
ring latch only 
center latch only 
both latches 
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Requiring center-ring coincidences with 

tron coincidence rate of 2/hr for both j “ Vt h e added qudifiedion 

signal of -0.6 ±0.4 The same comments apply as to Thble 1 , wnn 

-zzx. f jrr-r^xs 

- - - - 

lost when averaging over many runs or for a long duration. 

fte^t^expudnMtd 8 tonllops tt TWsaUOT^an^ai^““^ 1 ^^^' 

‘s^TtonneutonsingK snd which offer, .notor opportunity 

to 

system was incabable of accepting back-to- . —rents* the digitizing time for 

the electronics is specified* dedgned to “*^3 ps.K this was 
ADC’s and TDC’s will preclude ncquisriiou «f (<x ^tiple neutrons 

the case, then the only way to recogmze urs system was probably incapable , 

within an event. Since the event gdcvndthwss f et C „f(^ 

V. detecting an entire burst of longer ^ ^ ^ , t 

3 ^ Since one of the purposes of th xperim understood. I have written a 

is important that the response to multiple e vent ^ “ wious neutron burst 

MonteCarlo progmm to compute to proba cy totnmy Cotton. where 

multiplicities. The results have been checked for con _ ^ y ^ rountOT 

possible. The cdculdkms «n ’ * • fLdSe that to efficiencies 
of 0.8% efficiency as described in the Y proper expression, verified by 

predicted in the report i^ gr^sly ™ x neutron, 

where (“) is the usud binomid ^ssfn^S 

can make to expression much mm «mph < ^ opprojdnmtion to the above 

can become much greater than one, an j nV olve terms with Me, which can 

expression. Note that approximations for small e will involve terms 

not be assumed to be small. v e uniouely applied to a system 

In addition, the egression used m the report detectors could be 

with detectors of differing effidency. For examp t j iese reasons, the question 

for 10 ring and 1 central or 9 ring and 2 central the 

is most reliably handled by a Uontei advisee to limit 

r"irion“ detec toTmultipli city to tot for to ring counters. 
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„ /a\ _ n - 0.08) 70 = 18.9% 

P20,10,0.00#V U ; — V A 1 

Mid 125 neutrons ere listed in Appendix B. Th conditions on multiple hits m 

st?r» *-« *-fws: 

0 , 0 . 1 . ** 

‘Inefficiencies’ below). « not required’, ‘any’ hit, and ‘sngl hit 

The conditions for the central cou * d t()F discrimination. The ring coun ^ 

only. The last isfor cases more than one hit are excluded, and only , 

bein^ wwpt^° ThT ring-hit V^km N hits or g^t« 

In addition to probabilities for N hits, t e ghown at the bottom. 

£. *J*£ The edition The criteria «**. 

Equation (1) U -not required', end ring counter, ere only. 

Z R X" '^ =, 10 'eadiotring end cental det^. (Thu, the 

and the ‘any’ probabilities.) . . . „ advertized as 90ns, which is ex- 

The deadtime of the pulse-shape Wth ^ event gate of lO/is., this-woul 

tremely ^conaide^g th«r gc^i»^ “ tte ** gate, the inefficency 

■^rsrcsiS- - *■ tzzttzzt&Z 

TSlSXZ* «. — — * - - 

' detectors. 
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Multiplicity Data , vielded hundred events per hour 

Requiring latches for center-ring coin ^ en y A . • pu lse-height/pulse-shape cuts 
with ring multiplicities from o as g ‘ frQm ^ t h foreground and background, 

f " events mentioned in the .port ere duo to the looser 

CUtS Z S. of vuriouo applied cuts on the hourly data 

rate in ring detectors. 

Thble 2: CRNTM MM 

^4UUU 

~ 100-200 

~200 

~10 

~2 

0 . 5 ± 0.10 (28 evts) 


XX 

XX 

XXX 

xxxx 

XX 

XX 

XXX 

XXX 


2 

3 

2 

3 


0.25±0.08 (14 evts) > 


0 


where 

C is Center-detector latch 
R is Ring-detector latch 

I - rVt « *r\ 


R is Ring-detector iarcu 
N is Neutron pulse-shape/pulse height cut 
T is loose Time-of-fiight cut 


£> 


•r 


Bg(/hr) 

•>4000 

~ 100-200 

~200 

~10 

~2 

0.28x0.03 (7 evts) 

0 

0.1±0.05 (3 evts) 

0 


x lo WVWV# ****'' ~ Q 

M is Multiplicity for specified cut „ f for ,u good-lateh combinations 

Figure 8 . shows the distribution of tunes f g foreground M 4 events) and 

for all events with neutron multr pilot ws o^ne °y^“j ulicitieg ) points because some 
background (8 events). There are more than 17 (tore P ^ ^ time d; s tri- 

aims sx ksix" “ 

because only one start time from each cen groups.) This is a serious problem, 

that satisfied the CRN2 cut only 2 chance that the second 

because it means that at m um u*sm .. g rgt j ts ^tch will not be set, and time 
center counter will fire more than 500ns. aft » , ^ uo This is an example of 

information will not be avai} ab ^ e for ^Xlv^Xta without logic or electronics diagrams. 
‘ he "So^on^rthflter counters is aim not useful at high multiplicities 

^g^ 

%££££* * (ph - 720 - ps = 1300) - 

Figure 4. 
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the apparatus consists not only of the probability that a neutron willgen. 
crate appropriate pulses but that such p^s “oHh.‘ «£tro“re 

4V that lhe T er ; 

Wnt wi Lsembled and conducted in haste, there is some grounds for concern about 
rr“r iC s actually behaved as advertised. The following observatrons me 

pertmOTt. e l eclr0 nics diagram means that one can not be sure what correla- 

- io b“ expected and L proper in the data and waste nm^memfe™* 
thein. Some problems cannot be solved mid may not even be recogmzed without such 

. Aia^hthe logbook mentions that the detector gains were Rusted, the I«W«^ 
soectmshow great differences between detectors. Amplmer saturation was a Pattern. 

Thresholds, v£fch are of fundamental importance in any neat™ 

appear to^ry among the detectors. Figure 5 shows raw spectra from Cf runs which 

illustrate these problems. .. 

s , t W pre not adiusted, and several were at or near zero, making 

• makes it difficult to determine the eiectronic 

inefficient of the detector, which prefers an unambiguous electronic zero. 

• A small fraction of events had no latches (Fig. 7), which althoug pro y no 
significant is disquieting because such events should be logically impossible, and when 
SJJStfS is as low as in this experiment, any correlation between sudr losses 

and real data could be fatal. _ , 

• The TDC’s factually TAC’s) and at least some of the ADC’s were not cleared ef er 
each event Thus a detector without a latch set would contain data from a Previous 

!^t. This could be a problem if there areany spurious «>?£ 

■leftover’ events, which can be identified by comparing unlatched TDC values with 

the most recent latched value. . . ... , , • 

. There was substantial crosstalk between TDC’s, which was not entirely eliminated by 
requiring latches. Figures 2 and 8 also illustrates this. , 

• The TDC’s were not matched. This was actually helpful in identifying crosstalk, but 

demonstrates lack of care in setting up. 

• There was an enormous (about 30 times the nominal count rates) noise source m 
detector D1 which exactly simulated a good time of flight. Perhaps i was a p 
“merely or eccid it elly left in. In mtf c«e, it was only identrfiable by lack of 

. Tted^ln did not include any mechanism for determining the relative times of events 
in the two groups of detectors. Useful multiplicity data require latches forboth central 
del if Zv corresponding ring counters are included. This puts a 500ns window 
o:Tdo™3 P preba«y plains why so few op-down 
when center latches were required. A complete understan ng 
latches on multiplicity requires more MonteCarlo work. 
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Problems of this type may be expected from a new setup, but for a previously existing 
apparatus put together by asupposedly experienced group, they do not mcrease confidence 
that any results will be trustworthy. 

An obvious source of electronic inefficiency was that a substantial fraction of even s 
which had good times of flight failed the neutron-gamma discrimination cut 
had no pulse shape and/or no pulse height. Interestingly enough, for runs after run 30, 
when a|ustment of thresholds, etc. apparently ceased, the inferred ineffiomcy for da a 
mid background runs held sternly at about 60% (after neglecting the no.se in Dl), while 

‘ he ?™Sts h0 W “ 2SE.-S nries. the correction to the effiriencydue* 
findtime depends on what the correction was for the specific Cf runs used. Witho 
information on what actually went into the efficiency caclulations I cannot choosea 
rdtim:“on BaL onT observed variation both in the 
detection rates, I am indined to be suspicious that a single value of 0.8% is . aa 
apply to all. Detector U3, for instance seemed to have less than half the efficiency som 

° f th Thfrrrection for the more stringent neutron cut. that I used should probably be 

a b ° U ht I attribute the shove mentioned loss to deadtime in the pulse-height or 
aooaratuT then it is tempting to infer that the backgrounds have a somewhat lugher 
imfltiolicity (of gammas or muons) than that of the Californium (2-3 neutrons). 

The probability of self veto sms given as 0.6%, so the probability of veto for a burst 
of 100 neutrons is 45%, which could constitute a substantial inefficiency. 

The probability that a burst of 20 neutrons should have a rmg multiphaty of 3 or more, 
is from 5% to 12%, depending on the exact effect of deadtime. If afew such events oocurrah 
they wouM not be expected to be detected. Bursts of 125 neutrons have probabihties of 
from 78% to 99% of Mug detected. We can conclude that there were probably no events ^ , 

^°“je tw 50 neutrons in the data runs, although calculating the exact confidence l«ds J 
remains tobe done. The excess of dmible hits in the forced over he^poimd m 
Table 1, might be consistent with some smaller bursts, but it is difficult to set an upper 

limit on bursts of less than 50 neutrons. 

I have looked at a selection of the data consisting of about 400K f^FOund events 
representing about 56 hours worth of data. The selection also includes a corresp £ 
SXTound and mm. “’Cf (neutron source) data. No “Co (enmgy 
data wore available. This subset of the data contains the events mentioned m the mtenm 
report, and it is unlikely that subsequent review of the entire data set wih be 
change any of the conclusions. 


£yt /uJofa 
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pulse shape DO 



52BG 

Figure 1. Pulse shape vs. height in DO for background run 52 
. Distribution of events within neutron gate indicates they are 
not spill from gammas. 
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crosstalk 


no latches 
no center latch 
no ring latch 
both latches 


bad CfcR 
bed C, gammoR 

gamma C.gamme R 
leftove 

crosstalk 
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49 


gamma C 
bed R 



no latches 


no center latch 


no ring latch 
both latches 


tofz k pld#3 
52B6 

Figure 2. Particle ID vs time of flight for detectors D1 (#6) (© 
2.5ns/chan), and Do (10) (@256ns./chan) in background run 52. 
Part. ID is described in the text. 
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Figure 3. Distribution in time of center-ring 
coincidences as a function of the neutron mul- 
tiplicity, for all foreground and all background 
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location of 3rd ring neut. 
from Yale analysis 204Q 
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pulse shape U 

U#5 
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Figure 4. Pulse height spectrum for XJ5 in 
data run 45, showing location of the triple- 
neutron coincidence event which was rejected 
for having too large a pulse height 
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in 

34CF 

Figure 5. Pulse shape vs. height for all ring 
detectors for Californium run 34, showing 
mismatched gains and thresholds, amplifier 
saturation, and evidence of deadtime losses. 
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0-3 no letch: 
4-7 letch. 
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0-3 no 
4-7 lotch. 
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ph pld (sil) # 6 
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... 

(&). Data is from Californium run 34. 
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up center letch 
down center letch 
uo rl no letch mult 
down ring letch mult 

figure 7. Latch multiplicities for all foreground and background runs. 

Note events with no latches. 
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Figure 8. M scale (10 fts.) time of flight vs. detector 
for Californium run 34. Note nonuniform delays, and 
crosstalk from detector 10. 


✓ 




0 


S$> 


A 


14 


Appendix A. Neutron counts as a fraction of total clean events. 


Fire*- row after run name is total events with both ring and center latch set, with 
time "of flight between channels 0 and 20 (relative to gamma peak), and satisfying 
the 'ok' event gate on pulse shape/pulse height spectra. Second row is subset whic 
also satisfy the 'neutron' gate. Third row is ratio of previous numbers. 


tlow “0, thigh “ 20 

run* ul u2 u3 u4 u5 dl d2 d3 d4 d5 
22CF * 00029 00023 00022 00025 00000 00022 00025 00029 00022 00035 
(xx) 00001 00000 00000 00002 00000 00000 00001 00000 00001 00003 

0.034 0.000 0.000 0.080 0.000 0.000 0.040 0.000 0.045 0.086 
30BG* 00007 00002 00008 00002 00002 00006 00003 00012 00001 00004 

(XX) 00000 00000 00000 00000 00000 00000 00000 00001 00000 00000 

0 000 0.000 0.000 0.000 0.000 0.000 0.000 0.083 0.000 0.000 
31BG* 00009 00020 00012 00010 00019 00020 00014 00023 00006 00019 

(XX) 00000 00000 00000 00000 00000 00000 00000 00000 00001 00000 

0 000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.167 0.000 
32BG: 00312 00332 00428 00431 00337 00327 00431 00533 00350 00487 
(xxl 00000 00000 00000 00001 00000 00000 00004 00003 00004 00005 

0.000 0.000 0.000 0.002 0.000 0.000 0.009 0.006 0.011 0,010 
34CF: 09329 08602 10304 06314 11385 13230 11462 12334 08119 16076 
(CF) 01590 01510 00381 01731 02505 00948 01517 02523 5 2 2 4 Jo*, 

1 0.170 0.176 0.037 0.274 0.220 0.072 0.132 0.205 0.202 0.127 

40BG:00068 00091 00156 00026 00150 00356 00171 00233 00169 00493 
(BG) 00001 00000 00001 00000 00001 00002 00002 00003 00004 0C010 

0.015 0.000 0.006 0.000 0.007 0.006 0.012 0.013 0.024 0.020 
41BG*00001 00004 00009 00002 00006 00094 00011 00013 00004 00163 
(BG) 00001 00001 00000 00000 00001 00000 00000 00001 00000 00004 

1.000 0.250 0.000 0.000 0.167 0.000 0.000 0.077 0.000 0.025 
424A-00002 00001 00004 00001 00002 00105 00024 00004 00002 00243 
(FG) 00000 00000 00001 00001 00002 00000 C0002 00001 00000 00003 
* 0.000 0.000 0.250 1.000 1.000 0.000 0.083 0.250 0.000 0.012 

434A*00135 00130 00241 00090 00325 00442 00459 00228 00338 00684 
(FG>00002 00001 00000 00001 00002 00002 00002 00001 00009 00008 

0.015 0.008 0.000 0.011 0.006 0.005 0.004 0.004 0.027 0.012 
44#B:00075 00052 00148 00007 00220 00363 00307 00135 00129 00505 
(FG) 00000 00000 00001 00000 00003 00000 00005 00003 00001 00007 

0.000 0.000 0.007 0.000 0.014 0.000 0.016 0.022 0.008 0.014 
4540:00098 00086 00177 00058 00249 00346 00332 00419 00264 00538 
(FG) 00000 00000 00000 00000 00000 00002 00001 00002 00005 00007 

0 000 0.000 0.000 0.000 0.000 0.006 0.003 0.005 0.019 0.013 
464D-00067 00049 00129 00009 00204 00448 00367 00492 00215 00528 
(FG) 00002 00001 00000 00000 00003 00001 00000 00006 00001 00006 

0.030 0.020 0.000 0.000 0.015 0.002 0.000 0.012 0.005 0.011 
47#D: 00012 00007 00036 00002 00019 00202 00167 00146 00041 00250 
(FG) 00000 00000 00000 00000 00001 00001 00000 00001 00000 00001 

0.000 0.000 0.000 0.000 0.053 0.005 0.000 0.007 0.000 0.004 
52BG* 00029 00076 00007 00033 00107 00152 00152 00141 00095 00226 
(BG) 00000 00002 00000 00000 00000 00000 00001 00001 00001 00003 

0.000 0.026 0.000 0.000 0.000 0.000 0.007 0.007 0.011 0.013 
53#G* 00024 00065 00004 00008 00089 00151 00143 00118 00112 00206 
(FG) 00000 00001 00000 00000 00002 00000 00001 00001 00002 00002 

0.000 0.015 0.000 0.000 0.022 0.000 0.007 0.008 0.018 0.010 
average n fractio 


0?012 o!oil <h 003 0.008 0.011 0.003 0.007 0.010 0.017 0.013 


sum ave 
00232 
00008 

0.034 

00047 

00001 

0.021 

00152 

00001 

0.007 

03968 

00017 

0.004 

107155 

16390 

0.153 

01913 

00024 

0.013 

00307 

00008 

0.026 

00388 

00010 

0.026 

03072 

00028 

0.009 

01941 

00020 

0.010 

02567 

00017 

0.007 

02508 

00020 

0.008 

00882 

00004 

0.005 

01018 

00008 

0.008 

00920 

00009 

0.010 


type n 


all n/all 

FG: 00108 012278 0.009 +- 0.0008 
BG: 00040 003238 0.012 +- 0.0020 
CF: 16390 107155 0.153 +- 0.0012 
xx- 00748 015172 0.049 +- 0.0018 


/< 2 * 


k <jh 
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Appendix B. Probabilities for Ring-detector Multiplicities 


burst — 20 neutrons, count-10000 ei-0.008, eO — 

Center: ignore, any 

Ring: any sngl only any sngl 

Probability for exact ring multiplicity 
0:0.1838 0.1873 0.1838 ~ A 


only 


0.1596 0.1629 0.1596 
0.3165 0.3173 0.3132 
0.2502 0.2489 0.2461 
0.1131 0.1117 0.1103 
0.0320 0.0308 0.0306 
0.0058 0.0056 0.0056 
0.0003 0.0003 0.0003 
0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 

1ViViVVWW 0.0000 0.0000 0.0000 

accumulated probability for n or more ring hits 

q - ~ « a aaa aaei a onnc f\ OHHC O 


1:0.3593 0.3604 0.3558 
2:0.2839 0.2824 0.2793 
3:0.1298 0.1282 0.1267 
4:0.0365 0.0352 0.0350 
5:0.0061 0.0059 0.0059 
6:0.0005 0.0005 0.0005 
7 : 0.0001 0.0001 0.0001 
8 : 0.0000 0.0000 0.0000 
9 : 0.0000 0.0000 0.0000 
10 : 0.0000 0.0000 0.0000 


1.0000 1.0000 0.9871 
1:0.8162 0.8127 0.8033 
2:0.4569 0.4523 0.4475 
3:0.1730 0.1699 0.1682 
4:0.0432 0.0417 0.0415 
5:0.0067 0.0065 0.0065 
6:0.0006 0.0006 0.0006 
7 : 0.0001 0.0001 0.0001 
8 : 0.0000 0.0000 0.0000 
9 : 0.0000 0.0000 0.0000 
10 : 0.0000 0.0000 0.0000 
Expected Multiplicity: 
1.4967 1.4838 1.4869 


0.8776 0.8776 0.8658 
0.7180 0.7147 0.7062 
0.4015 0.3974 0.3930 
0.1513 0.1485 0.1469 
0.0382 0.0368 0.0366 
0.0062 0.0060 0.0060 
0.0004 0.0004 0.0004 
0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


0.10, 10 dets 
sngl 

any sngl 




0.1412 

0.2744 

0.2200 

0.1009 

0.0278 

0.0053 

0.0002 

0.0001 

0.0000 

0.0000 

0.0000 

0.7699 

0.6287 

0.3543 

0.1343 

0.0334 

0.0056 

0.0003 

0.0001 

0.0000 

0.0000 

0.0000 


0.1441 

0.2753 

0.2187 

0.0995 

0.0269 

0.0051 

0.0002 

0.0001 

0.0000 

0.0000 

0.0000 

0.7699 

0.6258 

0.3505 

0.1318 

0.0323 

0.0054 

0*0003 

0.0001 

0.0000 

0.0000 

0.0000 


ddtm— 0 .10 

only 

0.1412 
0.2715 
0.2162 
0.0984 
0.0267 
0.0051 
0.0002 
0.0001 
0.0000 
0.0000 
0.0000 

0.7594 
0.6182 
0.3467 
0.1305 
0.0321 
0.0054 
0.0003 
0.0001 
0.0000 
0.0000 
0.0000 


1.4992 1.4858 1.4890 1.5024 1.4888 1.4924 


predict 

(ign,only) 

0.1887 

0.3282 

0.2440 

0.1018 

0.0263 

0.0044 

0.00C5 

0.0000 

0.0000 

0.0000 

0.0000 

0.8939 

0.7052 

0.3771 

0.1331 

0.0313 

0.0049 

0.0005 

0.0000 

0.0000 

0.0000 

0.0000 

1.4007 


burst - 20 neutrons, count-10000 ei-0.005, eO 

Center: ignore, an y \ 

Ring: any sngl only any sngl only 

Probability for exact ring multiplicity 
0:0.3487 0.3712 0.3487 0.3043 0.3240 0.3043 

1:0.4017 0.3959 0.3794 0.3537 0.3489 0.3342 

2:0.1899 0.1794 0.1735 
3:0.0507 0.0458 0.0445 
4:0.0083 0.0071 0.0070 
5:0.0006 0.0005 0.0005 
6 : 0.0001 0.0001 0.0001 
7 : 0.0000 0.0000 0.0000 
8 : 0.0000 0.0000 0.0000 
9 : 0.0000 0.0000 0.0000 
10 : 0.0000 0.0000 0.0000 


0.10, 10 dets ddtm-1 
sngl 

any sngl only 


.00 


0.1669 0.1579 0.1523 
0.0442 0.0395 0.0383 
0.0074 0.0063 0.0062 
0.0006 0.0005 0.0005 
0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


JLV i ^ • VV V/V V • V w ~ . . , . , 

accumulated probability for n or more ring hits 
0:1.0000 1.0000 0.9537 * A onn ‘ > A ooeo 

1:0.6513 0.6288 0.6050 
2:0.2496 0.2329 0.2256 
3:0.0597 0.0535 0.0521 
4:0.0090 0.0077 0.0076 
5:0.0007 0.0006 0.0006 
6 : 0.0001 0.0001 0.0001 
7:0.0000 0.0000 0.0000 
8 : 0.0000 0.0000 0.0000 
9 : 0.0000 0.0000 0.0000 
10 : 0.0000 0.0000 0.0000 
Expected Multiplicity: 

0 . 9704 0 . 9236 0 . 9343 


0.0913 

0.1055 

0.0520 

0.0138 

0.0028 

0.0002 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.8772 

0.8772 

0.8359 

0.2656 

0.5729 

0.5532 

0.5316 

0.1743 

0.2192 

0.2043 

0.1974 

0.0688 

0.0523 

0.0464 

0.0451 

0.0168 

0.0081 

0.0069 

0.0068 

0.0030 

0.0007 

0.0006 

0.0006 

0.0002 

0.0001 

0.0001 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.9728 

0.9251 

0.9350 

0.9906 


0.0976 

0.1041 

0.0490 

0.0125 

0.0023 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.2656 

0.1680 

0.0639 

0.0149 

0.0024 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0913 

0.0993 

0.0472 

0.0119 

0.0022 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.2520 

0.1607 

0.0614 

0.0142 

0.0023 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


predict 
(ign, only) 
0.3585 
0.3774 
0.1698 
0.0429 
0.0067 
0.0007 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 

0.9560 

0.5975 

0.2202 

0.0503 

0.0074 

0.0007 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.9166 
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burst *= 20 neutrons, count-10000 
Center: ignore, any 

Ring:any sngl only 

Probability for exact rxng multiplicity 


ei-0.008, eO - 0 


only 


. 10 , 10 
sngl 
any 


<lets ddtm-1 . 00 
sngl only 


0:0.1838 

1:0.3593 

2:0.2839 

3:0.1298 

4:0.0365 

5:0.0061 

6:0.0005 

7:0.0001 

8 : 0.0000 

9:0.0000 

10 : 0.0000 


0.2159 0.1838 
0.3705 0.3286 
0.2673 0.2418 


0.1134 

0.0278 

0.0047 

0.0004 

0.0000 

0.0000 

0.0000 

0.0000 


0.1042 

0.0266 

0.0045 

0.0003 

0.0000 

0.0000 

0.0000 

0.0000 


0.1596 0.1879 0.1596 
0.3165 0.3257 0.2892 
0.2502 0.2361 0.2138 
0.1131 0.0988 0.0906 
0.0320 0.0245 0.0233 
0.0058 0.0043 0.0042 
0.0003 0.0003 0.0002 
0.0001 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


0.0461 0.0545 
0.0938 0.0972 


0.0756 

0.0332 

0.0090 

0.0017 

0.0001 

0.0001 

0.0000 


0.0713 

0.0279 

0.0074 

0.0011 

0.0002 

0.0000 

0.0000 


1 UIU.UUUU v/.vv/w v.wwww . .. a 

accumulated Probability for n or more ting hits 


1.0000 1.0000 0.8898 
0.8162 0.7841 0.7060 
0.4569 0.4136 0.3774 
0.1730 0.1463 0.1356 
0.0432 0.0329 0.0314 
0.0067 0.0051 0.0048 
0.0006 0.0004 0.0003 
0.0001 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
Expected Multiplicity: 

1 J4967 1.3824 1.4110 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0.8776 0.8776 0.7809 
0.7180 0.6897 0.6213 
0.4015 0.3640 0.3321 
0.1513 0.1279 0.1183 
0.0382 0.0291 0.0277 
0.0062 0.0046 0.0044 
0.0004 0.0003 0.0002 
0.0001 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


0.0000 0.0000 
0.0000 0.0000 

0.2596 0.2596 
0.2135 0.2051 
0.1197 0.1079 
0.0441 0.0366 
0.0109 0.0087 
0.0019 -0.0013 
0.0002 0.0002 
0.0001 0.0000 
0.0000 0.0000 




burst - 50 neutrons. 

Center: ignore, - 

Ring: any sngl only <!???; 

Probability for exact ring multiplicity 


0.0000 0.0000 
0.0000 0.0000 

1.4992 1.3851 1.4138 1.5039 1.3860 

count-10000 ei— 0.008, eO 
any 

only 


0.0461 

0.0855 

0.0640 

0.0262 

0.0070 

0.0011 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.2300 

0.1839 

0.0984 

0.0344 

0.0082 

0.0012 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 


predict 

(ign,only) 

0.1887 

0.3282 

0.2440 

0.1018 

0.0263 

0.0044 

0.0005 

0.0000 

0.0000 

0.0000 

0.0000 

0.8939 

0.7052 

0.3771 

0.1331 

0.0313 

0.0049 

0.0005 

0.0000 

0.0000 

0.0000 

0.0000 


1.4183 1.4007 


0.10, 10 
sngl 
any 


det s ddtci- 0 .00 


sngl only 


0:0.0167 

1:0.0788 

2:0.1986 

3:0.2673 

4:0.2367 

5:0.1356 

6:0.0522 

7:0.0119 

8 : 0.0020 

9:0.0002 

10:0.0000 


0.0167 

0.0788 

0.1986 

0.2673 

0.2367 

0.1356 


0.0167 

0.0788 

0.1986 

0.2673 

0.2367 

0.1356 


0.0522 0.0522 
0.0119 0.0119 
0.0020 0.0020 
0.0002 0.0002 
0.0000 0.0000 


0.0167 0.0167 
0.0785 0.0785 
0.1979 0.1979 
0.2663 0.2663 
0.2356 0.2356 
0.1348 0.1348 
0.0516 0.0516 
0.0118 0.0118 
0 .0020 0.0020 
0.0002 0.0002 
0.0000 0.0000 


0.0167 

0.0785 

0.1979 

0.2663 

0.2356 

0.1348 

0.0516 

0.0118 

0.0020 

0.0002 

0.0000 


0.0167 0.0167 
0.0785 0.0785 
0.1979 0.1979 
0.2663 0.2663 
0.2356 0.2356 
0.1348 0.1348 
0.0516 0.0516 
0.0118 0.0118 
0.0020 0.0020 


XU l U • uuyu v . v/v/v v v.vw , . , ^ 

accumulated probability for n or more rin, hits 


6 : 1.0000 1.0000 1.0000 
1:0.9833 0.9833 0.9833 
2:0.9045 0.9045 0.9045 
3:0.7059 0.7059 0.7059 
4:0.4386 0.4386 0.4386 
5:0.2019 0.2019 0.2019 
6:0.0663 0.0663 0.0663 
7:0.0141 0.0141 0.0141 
8 : 0.0022 0.0022 0.0022 
9 : 0.0002 0.0002 0.0002 
10 : 0.0000 0.0000 0.0000 
Expected Multiplicity: 
3.3170 3.3170 3.3170 


0.9954 0.9954 0.9954 
0.9787 0.9787 0.9787 
0.9002 0.9002 0.9002 
0.7023 0.7023 0.7023 
0.4360 0.4360 0.4360 
0.2004 0.2004 0.2004 
0.0656 0.0656 0.0656 
0.0140 0.0140 0.0140 
0.0022 0.0022 0.0022 
0.0002 0.0002 0.0002 
0.0000 0.0000 0.0000 


0.0002 

0.0000 

0.9954 

0.9787 


0.0002 

0.0000 

0.9954 

0.9787 


0.9002 0.9002 
0.7023 0.7023 


0.4360 

0.2004 

0.0656 

0.0140 

0.0022 

0.0002 

0.0000 


0.4360 

0.2004 

0.0656 

0.0140 

0.0022 

0.0002 

0.0000 


0.0167 

0.0785 

0.1979 

0.2663 

0.2356 

0.1348 

0.0516 

0.0118 

0.0020 

0.0002 

0.0000 

0.9954 

0.9787 

0.9002 

0.7023 

0.4360 

0.2004 

0.0656 

0.0140 

0.0022 

0.0002 

0.0000 


3.3148 3.3148 3.3148 3.3148 3.3148 3.3148 


predict 

<ign,only) 

0.0155 

0.0672 

0.1289 

0.1435 

0.1026 

0.0493 

0.0161 

0.0035 

0.0005 

0.0000 

0.0000 

0.5272 

0.5117 

0.4445 

0.3155 

0.1720 

0.0694 

0.0201 

0.0040 

0.0005 

0.0000 

0.0000 

2.9171 
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burst - 50 neutrons, count-10000 ei-0.005, eO - 

Center : ignore, an Y 

Ring:any sngl only any sngl only 
Probability for exact ring multiplicity 


0.10, 10 dets ddtm-0 
sngl 

any sngl only 


.10 


0:0.0779 0.0820 0.0779 
1:0.2298 0.2345 0.2259 
2:0.3063 0.3062 0.2978 
3:0.2277 0.2238 0.2190 
4 : 0:1141 0.1108 0.1094 
5:0]0344 0.0336 0.0330 
6:0^0078 0.0072 0.0071 
7:0,0018 0.0017 0.0017 
8 : 0]0002 0.0002 0.0002 
9 :0 , ; 0000 0.0000 0.0000 
10 : 0,0000 0.0000 0.0000 


'0.0772" 0.0813 0.0772 
0.2284 0.2331 0.2245 
0.3046 0.3044 0.2961 
0.2269 0.2231 0.2183 
0.1136 0.1103 0.1089 
0.0341 0.0333 0.0327 
0.0077 0.0071 0.0070 
0.0018 0.0017 0.0017 
0.0002 0.0002 0.0002 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


accumulated probability for n or more ring hits 


0:1.0000 1.0000 0.9720 
1:0.9221 0.9180 0.8941 
2:0.6923 0.6835 0.6682 
3:0.3860 0.3773 0.3704 
4:0.1583 0.1535 0.1514 
5:0.0442 0.0427 0.0420 
6:0.0098 0.0091 0.0090 
7:0.0020 0.0019 0.0019 
8:0.0002 0.0002 0.0002 
9 : 0.0000 0.0000 0.0000 
10 : 0.0000 0.0000 0.0000 
Expected Multiplicity: 
2.2149 2.1862 2.1988 


0.9945 0.9945 0.9666 
0.9173 0.9132 0.8894 
0.6889 0.6801 0.6649 
0.3843 0.3757 0.3688 
0.1574 0.1526 0.1505 
0.0438 0.0423 0.0416 
0.0097 0.0090 0.0089 
0.0020 0.0019 0.0019 
0.0002 0.0002 0.0002 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


0.0507 

0.1546 

0.2037 

0.1533 

0.0751 

0.0228 

0.0045 

0.0015 

0.0002 

0.0000 

0.0000 

0.6664 

0.6157 

0.4611 

0.2574 

0.1041 

0.0290 

0.0062 

0.0017 

0.0002 

0.0000 

0.0000 


0.0532 

0.1582 

0.2032 

0.1511 

0.0727 

0.0221 

0.0043 

0.0014 

0.0002 

0.0000 

0.0000 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


6664 

6132 

4550 

2518 

1007 

0280 

0059 

0016 

0002 

0000 

0000 


0.0507 

0.1523 

0.1977 

0.1477 

0.0717 

0.0218 

0.0042 

0.0014 

0.0002 

0.0000 

0.0000 

0.6477 

0.5970 

0.4447 

0.2470 

0.0993 

0.-0.276 

0.0058 

0.0016 

0.0002 

0.0000 

0.0000 


predict 

(ign,only) 

0.0769 

0.2025 

0.2350 

0.1583 

0.0685 

0.0199 

0.0039 

0.0005 

0.0000 

0.0000 

0.0000 

0.7657 

0.6887 

0.4862 

0.2512 

0.0929 

0.0244 

0.0045 

0.0006 

0.0000 

0.0000 

0.0000 


2.2158 2.1870 2.1997 2.2140 2.1855 2.1973 2.0226 




burst - 50 neutrons, count-10000 ei-0.008. 

Center: ignore, any 

Ring: any sngl only any sngl only 
Probability for exact ring multiplicity 

A. a1a« /■»! f\ Aino a nin A A1 an A QQ 


0.10, 10 dets ddtm-0 
sngl 

any sngl only 


,10 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0,0167 0.0193 0.0167 
0:0788 0.0845 0.0763 
0.1986 0.2083 0.1906 
0.2673 0.2668 0.2498 
0:2367 0.2334 0.2198 
0.1356 0.1286 0.1221 
0.0522 0.0464 0.0451 
0.0119 0.0106 0.0105 
0 j0020 0.0019 0.0019 
0.0002 0.0002 0.0002 
0.0000 0.0000 0.0000 


0.0167 0.0193 0.0167 
0.0785 0.0842 0.0760 
0.1979 0.2075 0.1899 
0.2663 0.2659 0.2489 
0.2356 0.2323 0.2187 
0.1348 0.1278 0.1213 
0.0516 0.0458 0.0445 
0.0118 0.0105 0.0104 
0.0020 0.0019 0.0019 
0.0002 0.0002 0.0002 
0.0000 0.0000 0.0000 


accumulated probability for n or more ring hits 


0 : 1.0000 

1:0.9833 

2:0.9045 

3:0.7059 

4:0.4386 

5:0.2019 


6 : 0 . 

7:0 

8:0 

9:0 

10:0 


0663 

0141 

0022 

0002 

0000 


1.0000 0.9330 
0.9807 0.9163 
0.8962 0.8400 
0.6879 0.6494 
0.4211 0.3996 
0.1877 0.1798 
0.0591 0.0577 
0.0127 0.0126 
0.0021 0.0021 
0.0002 0.0002 
0.0000 0.0000 


0.9954 0.9954 0.9285 
0.9787 0.9761 0.9118 
0.9002 0.8919 0.8358 
0.7023 0.6844 0.6459 
0.4360 0.4185 0.3970 
0.2004 0.1862 0.1783 
0.0656 0.0584 0.0570 
0.0140 0.0126 0.0125 
0.0022 0.0021 0.0021 
0.0002 0.0002 0.0002 
0.0000 0.0000 0.0000 


0.0111 

0.0523 

0.1311 

0.1807 

0.1608 

0.0910 

0.0339 

0.0083 

0.0011 

0.0001 

0.0000 

0.6704 

0.6593 

0.6070 

0.4759 

0.2952 

0.1344 

0.0434 

0.0095 

0.0012 

0.0001 

0.0000 


0.0125 

0.0563 

0.1379 

0.1799 

0.1591 

0.0862 

0.0303 

0.0070 

0.0011 

0.0001 

0.0000 

0.6704 

0.6579 

0.6016 

0.4637 

0.2838 

0.1247 

0.0385 

0.0082 

0.0012 

0.0001 

0.0000 


0.0111 

0.0510 

0.1259 

0.1687 

0.1499 

0.0818 

0.0290 

0.0070 

0.0011 

0.0001 

0.0000 

0.6256 

0.6145 

0.5635 

0.4376 

0.2689 

0.1190 

0.0372 

0.0082 

0.0012 

0.0001 

0.0000 


predict 

(ign,only) 

0.0155 

0.0672 

0.1289 

0.1435 

0.1026 

0.0493 

0.0161 

0.0035 

0.0005 

0.0000 

0.0000 

0.5272 

0.5117 

0.4445 

0.3155 

0.1720 

0.0694 

0.0201 

0.0040 

0.0005 

0.0000 

0.0000 


Expected Multiplicity 


3 . 


3170 


3.2477 3.2773 3.3148 3.2453 3.2747 3.3204 3.2513 3.2772 2.9171 
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burst ■= 50 neutrons, count-10000 
Center -.ignore, a n Y 

Ring: any sngl only any s ?gl 
Probability for exact ring multiplicity 


ei-C . 005, eO 


only 


0.10, 10 dets ddtm-1 
sngl 

any sngl only 


.00 


0:0.0779 0.1108 0.0779 
1:0.2298 0.2759 0.2013 
2:0.3063 0.3056 0.2356 
3:0.2277 0.1952 0.1572 
4:0.1141 0.0839 0.0716 
5:0.0344 0.0221 0.0194 
6:0.0078 0.0052 0.0045 
7:0.0018 0.0013 0.0011 
8 : 0.0002 0.0000 0.0000 
9 : 0.0000 0.0000 0.0000 
10 : 0,0000 0.0000 0.0000 


0.0772 0.1100 
0.2284 0.2744 
0.3046 0.3036 
0.2269 0.1947 


0.0772 
0.2000 
0.2341 
0.1569 

0.1136 0.0835 0.0713 
0.0341 0.0219 0.0192 
0.0077 0.0051 0.0044 
0.0018 0.0013 0.0011 
0.0002 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


.LUiVfWV/V V.VVV/V v.wwww , . . _ 

iced. ulated probability for « or more ring hits 


0:1.0000 1.0000 0.7686 
1:0.9221 0.8892 0.6907 
2:0.6923 0.6133 0.4894 
3:0.3860 0.3077 0.2538 
4:0.1583 0.1125 0.0966 
5:0.0442 0.0286 0.0250 
6:0.0098 0.0065 0.0056 
7:0.0020 0.0013 0.0011 
8 : 0.0002 0.0000 0.0000 
9 : 0.0000 0.0000 0.0000 
10 : 0.0000 0.0000 0.0000 
Expected Multiplicity: 
2.2149 1.9591 2.0325 


0.9945 0.9945 0.7642 
0.9173 0.8845 0.6870 
0.6889 0.6101 0.4870 
0.3843 0.3065 0.2529 
0.1574 0.1118 0.0960 
0.0438 0.0283 0.0247 
0.0097 0.0064 0.0055 
0.0020 0.0013 0.0011 
0.0002 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 


0.0023 

0.0075 

0.0100 

0.0067 

0.0034 

0.0007 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0306 

0.0283 

0.0208 

0.0108 

0.0041 

0.0007 

0.0000 

0.0000 

0.0000 

0.0000 

O.O0OO 


0.0033 

0.0094 

0.0095 

0.0058 

0.0023 

0.0003 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0306 

0.0273 

0.0179 

0.0084 

0.0026 

0.0003 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0023 

0.0066 

0.0076 

0.0044 

0.0020 

0.0003 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0232 

0.0209 

0.0143 

0.0067 

0.0023 

0.0003 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


predict 
(ign, only) 
0.0769 
0.2025 
0.2350 
0.1583 
0.0685 
0.0199 
0.0039 
0.0005 
0.0000 
0.0000 
0.0000 

0.7657 

0.6887 

0.4862 

0.2512 

0.0929 

0.0244 

0.0045 

0.0006 

0.0000 

0.0000 

0.0000 


eO - 


2.2158 1.9597 2.0338 

burst - 125 neutrons, count-10000 ei-0.008. 

Center: ignore, » n y 

Ring: any sngl only any sngl only 

Probability for exact ring multiplicity 

0 : 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

1*0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 

2:0.0042 0.0042 0.0042 0*0042 0.0042 0.0042 

3:0.0255 0.0255 0.0255 0.0255 0.0255 0.0255 


2.1144 1.8464 1.9181 2.0226 


0.10, 10 dets ddtiw-0 
sngl 

any sngl only 


4:0.0785 0.0785 0.0785 
5:0.1736 0.1736 0.1736 
6:0.2520 0.2520 0.2520 
7:0.2481 0.2481 0.2481 
8:0.1516 0.1516 0.1516 
9:0.0559 0.0559 0.0559 
10:0.0103 0.0103 0.0103 


0.0785 0.0785 0.0785 
0.1736 0.1736 0.1736 
0.2520 0.2520 0.2520 
0.2481 0.2481 0.2481 
0.1516 0.1516 0.1516 
0.0559 0.0559 0.0559 
0.0103 0.0103 0.0103 


accumulated probability for n or more ring hits 


0 : 1.0000 1.0000 1.0000 
1 : 1.0000 1.0000 1.0000 
2:0.9997 0.9997 0.9997 
3:0.9955 0.9955 0.9955 
4:0.9700 0.9700 0.9700 
5:0.8915 0.8915 0.8915 
6:0.7179 0.7179 0.7179 
7:0.4659 0.4659 0.4659 
8:0.2178 0.2178 0.2178 
9:0.0662 0.0662 0.0662 
10:0.0103 0.0103 0.0103 
Expected Multiplicity: 
6.3348 6.3348 6.3348 


1.0000 1.0000 1.0000 
1.0000 1.0000 1.0000 
0.9997 0.9997 0.9997 
0.9955 0.9955 0.9955 
0.9700 0.9700 0.9700 
0.8915 0.8915 0.8915 
0.7179 0.7179 C.7179 
0.4659 0.4659 0.4659 
0.2178 0.2178 0.2178 
0.0662 0.0662 0.0662 
0.0103 0.0103 0.0103 


0.0000 

0.0003 

0.0042 

0.0255 

0.0785 

0.1736 

0.2520 

0.2481 

0.1516 

0.0559 

0.0103 

1.0000 

1.0000 

0.9997 

0.9955 

0.9700 

0.8915 

0.7179 

0.4659 

0.2178 

0.0662 

0.0103 


0.0000 

0.0003 

0.0042 

0.0255 

0.0785 

0.1736 

0.2520 

0.2481 

0.1516 

0.0559 

0.0103 


0000 

,0000 

,9997 

,9955 

,9700 

.8915 

,7179 

.4659 

.2178 

.0662 

.0103 


0.0000 

0.0003 

0.0042 

0.0255 

0.0785 

0.1736 

0.2520 

0.2481 

0.1516 

0.0559 

0.0103 

1.0000 

1.0000 

0.9997 

0.9955 

0.9700 

0.8915 

0.7179 

0.4659 

0.2178 

0.0662 

0.0103 


.00 

predict 
{ign, only) 
0.0000 
0.0003 
0.0016 
0.0045 
0.0063 
0.0105 
0.0091 
0.0054 
0.0021 
0.0005 
0.0000 

0.0423 

0.0423 

0.0419 

0.0404 

0.0359 

0.0276 

0.0171 

0.0080 

0.0026 

0.0005 

0.0000 


6 . 3348 


6.3348 6.3348 6.3348 6.3348 6.3348 5.11-4 
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0 : 0.1 

1 : 0.1 


burst - 125 neutrons, 
de nter : ignore, 

Ring: any sngl onl 
~ '^ability for exact r 
I.Q000 0.0001 0.0000 
.0003 0.0013 0.0003 
2:0.0042 0.0097 0.0032 
0.0416 
0.1138 
0.2053 
0.2508 
0.2173 
0.1166 
9:0.<}559 0.0374 0.0341 
10:0.(103 0.0061 0.0061 


count-10000 ei-0.008, eO - 0.10, 10 dets ddtm-0.10 


3:0.0255 
4:0.0785 
5:0.1736 
6:0.2520 
7:0. S 481 
8:0.: 516 


0.0211 

0.0611 

0.1279 

0.1764 

0.1685 

0.0982 


any 

y any sngl only 

sngl 

any sngl 

only 

predict 

ing multiplicity 

0.0000 0.0001 0.0000 

0.0000 

0.0001 

0.0000 

(ign, onl 
0.0000 

0.0003 

0.0013 0.0003 

0.0001 

0.0004 

0.0001 

0.0003 

0.0042 

0.0097 0.0032 

0.0014 

0.0032 

0.0010 

0.0016 

0.0255 

0.0416 0.0211 

0.0074 

0.0126 

0.0061 

0.0045 

0.0785 

0.1138 0.0611 

0.0264 

0.0361 

0.0207 

0.0083 

0.1736 

0.2053 0.1279 

0.0543 

0.0640 

0.0400 

0.0105 

0.2520 

0.2508 0.1764 

0.0796 

0.0784 

0.0570 

0.0091 

0.2481 

0.2173 0.1685 

0.0765 

0.0699 

0.0538 

0.0054 

0.1516 

0.1166 0.0982 

0.0472 

0.0352 

0.0299 

0.0021 

0.0559 

0.0374 0.0341 

0.0191 

0.0130 

0.0122 

0.0005 

0.0103 

0.0061 0.0061 

0.0032 

0.0023 

0.0023 

0.0000 

’ for n or 
i 1.0000 

more ring hits 
1.0000 0.6969 

0.3152 

0.3152 

0.2231 

0.0423 

i 1.0000 

0.9999 0.6969 

0.3152 

0.3151 

0.2231 

0.0423 

; 0.9997 

0.9986 0.6966 

0.3151 

0.3147 

0.2230 

0.0419 

0.9955 

0.9889 0,6934 

0.3137 

0.3115 

0.2220 

0.0404 

1 0 .9700 

0.9473 0.6723 

0.3063 

0.2989 

0.2159 

0.0359 

0.8915 

0.8335 0.6112 

0.2799 

0.2628 

0.1952 

0.0276 

1 0.7179 

0.6282 0.4833 

0.2256 

0.1988 

0.1552 

0.0171 

> 0.4659 

0.3774 0.3069 

0.1460 

0.1204 

0.0982 

0.0080 

1 0.2178 

0.1601 0.1384 

0.0695 

0.0505 

0.0444 

0.0026 

: 0.0662 

0.0435 0.0402 

0.0223 

0.0153 

0.0145 

0.0005 

0.0103 

0.0061 0.0061 

0.0032 

0.0023 

0.0023 

0.0000 

! 6.3348 

5.9835 6.2352 

6.3350 

5.9971 

6.2474 

5.1144 


Oil.dOOO 1.0000 0.6969 
1:1.(000 0.9999 0.6969 
2:0. <997 0.9986 0.6966 
3:0. <955 0.9889 0.6934 
4:0.^700 0.9473 0.6723 
5:0.8915 0.8335 0.6112 
6:0.7179 0.6282 0.4833 
7:0.4659 0.3774 0.3069 
8:0.2178 0.1601 0.1384 
9:0.0662 0.0435 0.0402 
10:0.0103 0.0061 0.0061 
Expected Multiplicity: 
6.3348 5.9835 6.2352 


burst — - 125 neutrons, count— 10000 el— 0.008, eO 
Center : ignore, any 

Ring: any sngl only any sngl only 
Probability for exact ring multiplicity 



0.10, 10 dets ddtm-1 
sngl 

any sngl only 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0.0000 0.0088 0.0000 
0.0003 0.0577 0.0003 
0.0042 0.1548 0.0017 
0.0255 0.2442 0.0034 
0.0785 0.2417 0.0074 
0.1736 0.1701 0.0108 
0.2520 0.0894 0.0090 
0.2481 0.0271 <$.0042 
0.1516 0.0053 0.0016 
0.0559 0.0009 0.0004 
0.0103 0.0000 0.0000 


0.0000 0.0088 0.0000 0.0000 0.0000 0.0000 


0:1.0000 1.0000 0.0388 
1:1.0000 0.9912 0.0388 
2:0.9997 0.9335 0.0385 
3:0.9955 0.7787 0.0368 
4:0.9700 0.5345 0.0334 
5:0.8915 0.2928 0.0260 
6 : 0 .i7 1 7 9 0.1227 0.0152 
7:0.M659 0.0333 0.0062 
8:0.2178 0.0062 0.0020 
9:0.0662 0.0009 0.0004 
10:0.0103 0.0000 0.0000 
Expected Multiplicity: 
6^3348 3.6938 5.0851 


0.0003 

0.0577 0.0003 

0 .0000 

0.0000 

0.0000 

0.0042 

0.1548 0.0017 

0.0000 

0.0000 

0.0000 

0.0255 

0.2442 0.0034 

0.0000 

0.0000 

0.0000 

0.0785 

C.2417 0.0074 

0.0000 

0.0001 

0.0000 

0.1736 

0.1701 0.0108 

0.0000 

0.0000 

0.0000 

0.2520 

0.0894 0.0090 

0.0000 

0.0000 

0.0000 

0.2481 

0.0271 0.0042 

0.0000 

0.0000 

0.0000 

0.1516 

0.0053 0.0016 

0.0001 

0.0000 

0.0000 

0.0559 

0.0009 0.0004 

0.0000 

0.0000 

0.0000 

0.0103 

0.0000 0.0000 

0.0000 

0.0000 

0.0000 

ir n or 

more ring hits 




1.0000 

1.0000 0.0388 

0.0001 

0.0001 

0.0000 

1.0000 

0.9912 0.0388 

0.0001 

0.0001 

0.0000 

0.9997 

0.9335 0.0385 

0.0001 

0.0001 

0.0000 

0.9955 

0.7787 0.0368 

0.0001 

0.0001 

0.0000 

0.9700 

0.5345 0.0334 

0.0001 

0.0001 

0.0000 

0.8915 

0.2928 0.0260 

0.0001 

0.0000 

0.0000 

0.7179 

0.1227 0.0152 

0.0001 

0.0000 

0.0000 

0.4659 

0.0333 0.0062 

0.0001 

0.0000 

0.0000 

0.2178 

0.0062 0.0020 

0.0001 

0.0000 

0.0000 

0.0662 

0.0009 0.0004 

0.0000 

0.0000 

0.0000 

0.0103 

0.0000 0.0000 

0.0000 

0.0000 

0.0000 

6.3348 

3.6938 5.0851 

8 .0C00 

4.0000 

0.0000 


.00 

predict 

(ign,only) 

0.0000 

0.0003 

0.0016 

0.0045 

0.0083 

0.0105 

0.0091 

0.0054 

0.0021 

0.0005 

0.0000 

0.0423 

0.0423 

0.0419 

0.0404 

0.0359 

0.0276 

0.0171 

0.0080 

0.0026 

0.0005 

0.0000 

5. 1144 
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Appendix C. Inferred Electronic Efficiencies. 


Denominator (sum) is the sum of those events with a latch for the given detector and 
its corresponding center detector and time of flight between 0 and 20 channels, 
(relative to gamma peak). Numerator (ok) is subset which also satisfied the 'ok' ev 
gate on pulse shape/pulse height spectra. D1 is not included in the averages becaus 
of the presence of obviously spurious events in its denominator. 

run : ul 

22CF: 0 . 66 
23?? : 0 . 87 
30BG: 0 . 89 
31BG: 0.(57 
32BG: 0 .72 
34CF:0.86 
35PB:0.47 
36PB:0.25 
40BG:0.33 
4 1BG : 0 . 30 
42#A:0.35 
43#A:0.30 
44#B:0.29 
45#C:0.28 
46#D:0.25 
47#D:0.23 
52BG:0.28 
53#G:0.24 


u2 

u3 

u4 

u5 

dl 

d2 

d3 

d4 

d5 


ave*. 

sum*, ok* (*w/o 

0 • 80 

0.42 

0.56 

0.00 

0.49 

0.56 

0.56 

0.54 

0.59 

-> 

0.521091 

000432 

000247 

0.95 

0.68 

0.75 

0.97 

0.62 

0.72 

0.74 

0.85 

0.76 

-> 

0.809598 

014939 

011685 

1.00 

0.64 

0.38 

1.00 

0.17 

0.50 

1.00 

1.00 

0.83 

-> 

0.804453 

000067 

000051 

0.81 

0.75 

0.58 

0.95 

0.08 

0.41 

0.55 

0.87 

0.48 

-> 

0.673284 

000257 

000160 

0.87 

0.66 

0.72 

0.93 

0.08 

0.47 

0.53 

0.83 

0.55 

-> 

0.697263 

006159 

003997 

0.94 

0.68 

0.77 

0.96 

0.63 

0.73 

0.74 

0.86 

0.77 

-> 

0.812243 

184127 

147378 

0.70 

0.65 

0.64 

0.89 

0.04 

0.34 

0.35 

0.87 

0.46 

-> 

0.596772 

001154 

000581 

0.65 

0.64 

0.40 

0.97 

0.03 

0.35 

0.33 

0.81 

0.48 

-> 

0.541490 

000401 

000199 

0.53 

0.62 

0.44 

0.86 

0.04 

0.27 

0.34 

0.85 

0.46 

-> 

0.522281 

006604 

003010 

0.53 

0.62 

0.61 

0.86 

0.03 

0.31 

0.34 

0.86 

0.44 

-> 

0.540300 

002681 

001236 

0.53 

0.64 

0.59 

0.88 

0.03 

0.33 

0.34 

0.82 

0.40 

-> 

0.540387 

004750 

002167 

0.53 

0.61 

0.58 

0.88 

0.03 

0.30 

0.36 

0.84 

0.44 

-> 

0.538271 

009573 

004418 

0.56 

0.62 

0.43 

0.89 

0.03 

0.31 

0.36 

0.84 

0.45 

-> 

0.527964 

006969 

003244 

0.54 

0.62 

0.63 

0.90 

0.04 

0.28 

0.38 

0.83 

0.43 

-> 

0.541859 

007110 

003256 

0.47 

0.62 

0.48 

0.87 

0.03 

0.25 

0.33 

0.80 

0.45 -> 

0.503666 

007735 

003294 

0.51 

0.61 

0.58 

0.87 

0.03 

0.30 

0.34 

0.85 

0.45 -> 

0.527601 

003945 

001793 

0.60 

0.59 

0.62 

0.86 

0.04 

0.29 

0.29 

0.85 

0.46 

-> 

0.536787 

002752 

001269 

0.52 

0.61 

0.51 

0.89 

0.03 

0.29 

0.33 

0.88 

0.45 

-> 

0.525211 

002639 

001237 
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To: Yale/BNL/BYU collaboration 

From: S. E. Jones 

Re: Analysis of Yale data from August run 

Date: 26 December 1989/revised 3 January 1990 


Dear Colleagues: 


We set out in the collaboration to search for neutron bursts 
of the type reported by Menlove et al . of the LANL/BYU 
collaboration, and specifically to determine the time structure of 
neutron bursts if any were found. Since receiving the data tapes 
on October 26, the day before the Interim Report (IR) was submitted 
to the DOE/ERAB committee on cold fusion, we have found in the data 
and in the IR a number of items to report to the collaboration. 
The analysis was done with the help of Stuart Taylor, Russell 
Hunter and Alan Anderson. Al has worked with me on analysis of 
neutron data from muon catalysis experiments for years, including 
interpreting information from liquid-scintillator counters as also 
used at Yale, and has written up a separate report which is 
attached. Notice that the term "TDC" is used for "time to digital 
conversion" though at Yale the TDC was built up from TAC's rather 
than with the use of TDC units per se. 

1. The two central counters have the highest efficiency (reported 
in IR at approx. 10% each), yet these counters were not timed 
relative to each other (as I learned in October from Kelvin). Why 
not? These counters give us the best chance of detecting neutrons 
and since we are after time information, the lack of this 
information for the most efficient counters should be explained. 
This oversight (?) is not fatal IF we have good efficiency for 
burst detection in the ring counters AND reliable time (TDC). 
information between the central and ring counters. Unfortunately, 
these conditions are not well met in the Yale set-up, as we shall 
see. We are not happy with this situation; after analyzing it, we 
propose a solution based on what we have learned. 

2. A major problem in the current experiment is the low efficiency 
for neutron-burst detection coupled with the short running time 
(103 hours of foreground data, 28 hours of background). Al 
demonstrates in his report that the efficiency of the set-up for 
detecting neutron bursts is much less than reported in the Interim 
Report. The burst-detection efficiency seemed suspiciously high 
to me and the efficiency for detecting two neutrons in one ring 
counter suspiciously low. Al flagged the "100%" efficiency 
statements in the Interim Report: 


The most probable neutron multiplicity (Mq) can be calculated 
from an event of fold-(K) using the efficiency for detecting an event of 
fold K given by: 


«(K) - Hrf<c(12)*(M n -l)*<(ll)*...*(M n -K)* <(12-K) 

For neutron source of multiplicity 125 ve obtain 100% efficiency for 
detecting a neutron in each detector (fold-12), and for multiplicity of 20 

ve have 100% efficiency of detecting 5 neutrons. Hence a 5 (6) fold event 

most probably corresponds to a source with neutron multiplicity of M n - 20 

(• 24 ) . 


I 


The equation used m the Interim Report (above) t ttiit to • f *Ci€S 

trcf trebiro* th#* p tjv. 

Calculating the burst-detection efficiency can be done using 
Monte Carlo techniques to generate neutrons and to evaluate 
detection probabilities in the counters. A1 took the time to write 
such a Monte Carlo and proceeded to calculate the efficiency of the 
Yale set-up as a function of the multiplicity (M) of the burst. 

He also provides an analytical expression which he used to check 
the Monte Carlo. The results are given in tables in his write up. 

Note that the probability for detecting a burst of 125 neutrons in 
all ring counters is less than 5%, rather than 100% as stated 
erroneously in the Interim Report. Similarly, the probability of 
detecting 20 neutrons in 5 counters (2 central and 3 ring) is only 
about 15% (-6% when reasonable electronic inefficiency is folded 
in) — not 100% as stated in the Interim Report. 

On the other hand, the Monte Carlo calculation shows that the 
probability for detecting two neutrons in a ring counter for a 
large burst is not negligible (stated as 0.2% in IR). A related 
problem with low efficiency for burst-detection is that 
interpreting the initial multiplicity from the number of detectors 
fired is ill-defined, as shown in Al's tabulated probabilities. 

Thus, for example, the multiplicity 5 event in run 45 in the IR, 
could have originated from a burst of 50-100 neutrons, and cannot 
be rejected as a burst candidate just because counter U5 shows 
greater than 2.5 MeV energy. A burst this large also has a 
significant chance (-50% for 100 neutrons) of self -veto since there 
was no shielding placed between the veto counters and the sample 
cylinders. In particular, the Monte Carlo calculation shows that 
a burst of 50 neutrons has -25% probability of producing the 
counter-hit pattern of "event 45" including a double hit in U5. 

(See Table 1.) This is not negligible. It is worth noting that 
this event occurs at 90 minutes into warm-up from LN, temperature, 
corresponding to a sample temperature of about — 10°C. On the other 
hand, event 45 could easily represent a cosmic-ray spallation event 
of, sav 20 neutrons. We conclude that the low probabilities for 
burst detection coupled with low statistics (and other problems 
discussed by Al) preclude unambiguous evaluations of such events. 

With more data on time-structure of bursts, one might be able to 
distinguish cosmic-ray events from "real" bursts. As we discussed 
months ago, cosmic-ray induced particles JJ9#$f1v appear within 
nanoseconds whereas any "real" events could well produce neutrons 
spread out in time , perhaps over microseconds . ( See Table 1 . ) 
(Problems with timing information are discussed in parts 1 and 3.) 

One might be able to compensate somewhat for the low burst- 
detection efficiency of the Yale set-up with ample statistics from 
very long running, but we certainly don't have it. Samples used 
at Yale were taken to Brookhaven for longer examination as we 
agreed; I will leave it to Kelvin to discuss results. We must not 
ignore these results . However , the BNL detector cannot give us the 
information we seek regarding the initial time distribution of 
neutrons in a burst. 
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If our analysis is correct, then the burst-detection 
efficiency given in the Interim Report is greatly overoptimistic 
and the conclusions that follow from it are misleading. In 
particular, the stated upper limits on neutron bursts are way too ;?c ,„/. t 
low. Let's look once again at the number of bursts of 20 or more/f ’ 
expected in the Yale set-up, based on the Henlove results and 
efficiency numbers from Al's report: 

Number of bursts expected at Yale: 


We multiply the time (T) at each lab by the mass (M) of 
material per measurement hour, (~200g at Yale, -lOOg at LANL) . The 
number of bursts detected (“LANL) is given in the Menlove Update 
sent out the collaboration and is 39 bursts of > 20 neutrons over 
13,000 measurement hours (at ~lOOg/hr)0 One may quibble over 

corrections to this calculation for instance, using time spent 

in the -80° to 0°C range instead of total foreground running time, 
and the number of bursts detected at LANL in this range, increases 
the expected yield by about 3. 

But the overriding factor is the ratio of efficiencies (E) for 
burst detection (“Yale/*LANL) . Our analysis shows that the Interim 
Report was misleading on this crucial point. In particular, “Yale 
drops off sharply at N < 100. For example, the efficiency for 
detecting 20 neutrons is -100% in the Menlove counter^ while only 
about 12% in the Yale set-up if we require 3 or more ring counters 
to be hit (not 100% efficient as stated in the Interim Report) . 
This drops to about 6% when reasonable inefficiencies in the 
electronics are folded in (see Al's report). Thus, the ratio of 
burst-detection efficiencies for > 20 neutrons in a burst is 
significantly less than 1, so that: 


3. In addition to the problem of low efficiency for neutron-burst 
detection in the Yale set-up, other difficulties should be reported 
because they show what can be improved in future experiments. 
Figs. 1 and 2 show what appear to be electronic artifacts/ 
saturation in detectors. While not necessarily fatal, these 
artifacts should be cleaned up or at least understood. Fig. 3 
shows the existence of considerable cross-talk between timing 
channels. It is evident that a stop in one counter may generate 
stops in other timing channels, which may invalidate the time 
information in these other units. Now, imposing latches cleans up 


"Yale = "LANL x “Yale x “Yale * “Yale 

“LANL X "LANL “LANL 


= 39 x 103 hrs x -20Qq each hr „ “Yale 
13,000 hrs x -lOOg each hr “LANL 






X, v' r ns ~ 

"Yale < 1 burst’expected, primarily because of the 
low efficiency for burst detection. 



1 1 Qn fru-rffc' <vf 
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the problem appreciably as shown in Fig. 3b. But this imposes the 
constraint that neutrons be DETECTED within a 500 ns window and 
reduces what we can say about bursts — we don't know the time 
structure of bursts in advance and should not bias the observations 
by imposing a short window on detection. I was assured in August 
that we could analyze the data without imposing the latch 
constraint, but this looks impossible with regard to timing 
information. Furthermore, there is a significant probability of 
multiple hits in counters for a burst event and there is ambiguity 
in the timing information between neutrons emanating from the 
target versus neutrons recoiling from other detectors . 


What can be done? Certainly the problems in the set-up and 
the Interim Report outlined here and in Anderson's report cannot 
be ignored and the IR submitted for publication with minor changes. 
We need to take what we have learned and fulfill our commitment to 
do a careful, definitive experiment. nXS* 

I therefore propose an improved follow-up experiment. I am 
willing to run again at Wright Laboratory if the improvements 
discussed herein are made and if sufficient running time is 
available. Otherwise, I propose an experiment using the muon- 
catalysis neutron detector at LAMPF. I have been a spokesman for 
muon-catalyzed fusion experiments at LAMPF for nearly eight years; 
we may use the system until May 1990 when the 800 MeV proton beam 
recommences. The LAMPF apparatus has been used successfully to 
detect 150 neutrons (average) produced over several microseconds 
by muon-catalyzed fusion and can be reconfigured for cold fusion. 
The electronics are already in place and can be used with minor 
modifications. Quantities of shielding materials are also 
available while the proton beam is off (until May). The time- 
structure of neutron events from muon catalysis has proven 
extremely informative and we have the tools needed to evaluate the 
time spectra (Phys. Rev. Lett. 56# 588 (1986), PRL 51, 1757 (1983), 
Nature 221, 127 (1986), etc.). 

The apparatus uses 5 liquid-organic scintillators in 5" diam. 
x 5" deep cylinders. We would probably configure 4 in a 
tetrahedral arrangement to get the best burst-detection efficiency 
with the fifth detector to detect backgrounds simultaneously. Of 
course, we would perform Monte-Carlo calculations (as we have done 
for our muon-catalysis work) before settling on the configuration 
of detectors/samples . The Monte Carlo would include effects of 
neutron recoil in the counters and geometrical effects stemming 
from the fact that the sample is extended rather than point-like. 
(Note that such a study has not been done for the Yale set-up.) 
A brief write-up on the electronics by A. J. Caffrey is enclosed 
to provide further information. 

Here are improvements expected relative to the experiment at 

Yale: 

1 — Good efficiency for detecting bursts of neutrons; the four or 


b 


five central counters will be timed together (this is already in 
the LAMPF electronics) . We will use a Monte Carlo to assess 
whether shielding is needed between counters (or higher energy 
thresholds) to permit differentiation of neutrons coming from the 
target versus neutrons recoiling from other counters. Shielding 
will be placed between the cosmic-ray veto and the sample cylinder 
to reduce the probability that burst neutrons (if any) will 
generate a false veto. 

2 — no cross-talk between TDC's. Gains between the 4 detectors 
are rather closely matched at LAMPF facilitating determinations of 
the efficiencies of the counters when various n/g cuts are later 
imposed. The n/g discrimination at low-energies at Yale is good, 
but we do well at LAMPF (see attached report by A. J. Caffrey) and 
the system is well understood and clean after seven years of use. 
The logic-router circuitry at LAMPF (see Caffrey report) improves 
handling of multiple hits in a single counter. In particular, we 
can handle up to eight neutrons per detector, using 8 separate 
ADC's and 8 TDC's, so that we circumvent the -80 jusec deadtime in 
ADC's at Yale. ADC's for pulse-height information are charge- 
integrating rather than peak-sensing. An electronics diagram will 
be provided. 

3 — we have approval to use the equipment until the end of April 
(the proton beam restarts in May 1990). We should plan at least 
2 weeks for the first set of experiments — say in Feb. 1990. We 
will want first to characterize the background using our on-line 
analysis package. Following data analysis, we can run again (in 
April) if necessary. Thus, we will have a longer time than at Yale 
in August and opportunity for follow-up experiments. The beam is 


No, it won't be easy or "quick and dirty." But by pooling our 
resources and building on what we have learned, we can together do 
the definitive experiment we set out to do. 

p.S. In addition to our current collaboration, I have invited 
Menlove — who was invited before but could not come — and experts 
on the LAMPF system (Anderson, Caffrey and Paciotti) who have all 
said they would be glad to participate should we decide to run at 
LAMPF. 

P.P.S. We may also be able to do experiments with deuterium plus 
TRITIUM. The needed safety equipment is in place at LAMPF from 
muon catalysis experiments. We have had years of experience with 
d-t fusion experiments and will need the help of Richard Maltrud 
of LANL for experiments involving tritium. 
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(A copy to Gaj ewski, as I recently learned from him that he has been getting 
copies of Jones commuincations to me in the past) . 


Dear Steve, 


V0- 


Thanks for sending to me a copy of your analysis of the Yale data. 

Unfortunately, this package does not represent an analysis of the Yale 
data, since I don't see any plots of bursts detections, random emission etc. 

. per run. Instead this package represent, as I shall detail below a false use 
of the Yale setup. I only wish you or your student would have called me to 
discuss it before you waisted such a long time. 

In addition unfortunately, the report prepared by your student is infested with 
misundersatnding of what was done, as I shall detail below. Again I wish he 
would have returned my calls and discussed with me his work, it would have 
saved a great deal of time for him and all of us. I understand that you are 
eager to show that the Yale experiment is flawed, but you should also realize 
that understanding the set up is crucial for analyzing the data. 

I shall reply point by point to your questions : 

1. You suggest that there is a flaw in the setup, since the time between the 
two central detectors is not measured. Let me then remind you that in order to 
study burst emission over 20 microseconds with a few nanosecods accuracy (as we 
set out to do) we were obliged to use separate TACs in conjunction with" ADCs 
for each detector . No existing TDC even at the high energy- electronics pool in 
Brookhaven would have allowed us to perform this task. Since we have in our 
laboratory only 6 TACs, a request was sent to all laboratories around us to 
provide us with the necessary electronics. You will recall that we could only 
get extra electronics from Stony Brook, Brookhaven and MIT-Bated linear 
accelerator. You promissed to look into bringing more TACs, but found none. 
/Since we could come up with only 12 TACs a decession was maid to have two 
Ajseparate hemispheres each working independently and to require that the timing 
y of the two central detectors are matched. foj ( J et ; s -ic ^ 
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3. The cross talk as I metioned to you is a well known problem of the Phillips 
Scientific discriminators. A similar effect is seen in many experiments at the 
Brookhaven AGS (e . g. E814) . The problem is well familiar and standard 
solutions exist. Unfortunately you chosed the worst solution, i.e. to gate on 
the latch. A better solution is to gate on information included in an ADC 
channel. Eor example if you have a completely open gate on the energy, i.e. 
all erier gies _accgpt£d__but you require a Pulse height information of any kind 
therj gio cross talk ' ex ist^X no data is lost and no stringent timing are imposed 
as in youF~solutIbfflSHIch\uses the latch. . _ . / . r 

^ Cross- 7c, I !< std J Sezfit'f fesfc*jr 

You also mentioned some problem with saturation. Again you should realize that 
in order to have as low threshold (of 50-70 keV electron equivalen) as we used 
in our experiment (in fact the lowest ever used on a neutron detector of that 
kind in an active research!), this requires driving the gain very high, as we 
did, which required saturating the channel at some 4.5 MeV or so. What is the 
problem with saturation at 4.5 MeV, if one looks at 2.5 MeV neutrons? 


2. Unfortunately your student used the wrong equation for hit probability. 

I am surprised you allowed his report to leave BYU in the form it did. So 
here are the comments on your student work: 

A. The equation your student used to calculate hit probability is wrong. 
Consider the following situation with one detector of 100% efficiency: for a 
neutron source of multilplicity 1 his formula gives hit probability of 100%, 
correctly. Eor multiplicity of 2 it yield a hit probability of 0%. This 
happens since he requires one and only one hit- and in the second case we have a 
100% chance of a double hit. If you inspect in a similar way his formula for a 
detector with 50% efficiency and again only one edetctor and multiplicity 
ranging from 1 to N, you again will notice that in fact as the mutiplicity 
rises the efficiency decraeses, as one may expect for requiring one and only 
one hit in a detector . This is certainly not the case in our experiment since 
we have no idea whether in fact the hits are single or multiple. In fact for 
the central detectors we are for sure seeing doubles etc . even with 
multiplicity of 20. Your student should read the paper that I mentioned to you 
over Bitnet last Friday and carefully correct this mistake. Correction of this 
mistake would be done by summing over all single and all multiple hits. 

Another problem in your student analysis is that he assumes (may be without 
even knowing) that , all detectors have the same efficiency. This is why in 
our report we don't give the formula with powers of epsilon but as e(i)s as the 
efficiency is i dependent (if all e(i)=e one retreave the exponent used by your 
student) . The math for this is somewhat involved but already published, for 
example in the paper I e-mailed to you. 

In the example we have above with one detector with efficiency of 100% and a 
neutron source of multilicity 2, using my equation yield an efficiency of 200%. 
In this sense it is simply interpreted as 100% probability for a double hit. 

In fact you will notice that the formula we wrote in our report corresponds to 
the first binomial times K!*e**K, and the rest of the factors in (1) need to be/ 
rewritten and summed over. While having efficiency larger than 100% as 
suggested in my formula are indeed strange, they simply have to be 
explained (i.e. they are not wrong) . 
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May I suggest that we discuss it 
trivia . 


;o as to avoid all this waist of time on 


In the report your student talk about some advertisement of veto efficiency. 
The veto efficiency was not advertized, it was measured! if your student need 
to understand this measurement I shall be happy to explain it, as detailed m 
the run book. If he is missing an electronic diagram he should get in couch 

with me and I will be hapr-v n : provide him with one. It is not clear at ai l 
from his reoort where w is is analysis hampered by not having the electronics 
diagram. Could you cinri . / i. r :.: 5 - statement? 



Again I am very eager to bring our data to publication and I am not sure why it 
cook ‘some two monts to get to the point where we discuss how to calculate the 
efficiency of our detector? I am very interested in releaving all problems of 
communications, I am surprised that my calls to your students where never 
returned and that the data analysis, is dragging so slowly, even though we are 
so far only discussing trivia. 

Best Regards Moshe gai. 
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Dear Steve, 

I have been trying very hard to get in touch with you to discusse the trivial 
mistakes that your student (A1 Anderson) makes in his attempt to analyze the 
Yale data. Steve the reprot you mailed to us is so bad that Pauli would have 
said: "it is hot even wrong!". You could have at least asked the student to 
use a more polite language! To hear from him that the Yale experiment was 
performed by amatures, after all what this laboratory did to allow this 
research is just not acceptable. I understand that you are now circulating 
this false criticism to whoever would listen to you. We learned about it from 
Kevin Wolf at Texas A&M from Gajewski at DOE etc. In fact Kevin Wolf knew 
about it before myself. Let me assure you that Yale University will not take 

part in this discussion! T hunst no* diu u i ^ Wolf - not _u ags 

ft ^ ^ ‘TH'iM* Oclob*/* L f€t cl \ & if)-/- 

Concerning the data analysis by your Student. I feel that it is so bad that Vc /. 

we might as well just ignore it. I will pretend that I never got it, in order 
to save you some embarrasment . However I would be happy to discuss it with 
you on the phone if you so desire. 

Since your student managed to confuse you on the issue of calculating hit 
probabilities, I would like to take the time to correct the algebra. I send to 
you by FAX two pages of algebra. You can see that once one corrects the wrong 
formula used by your student and with a little bit of algebra one gets excatly 
the formula of our report times a factor: (M-K) /K! ^sThis factor is of the order 
unity in our case (M=12 andK=3) so it was neglected \in the report. We can add 
it in the final publication—^ ^ rt- 10 *** *"'*’*? 

9 n- 7oO to nl OV C /(tju 
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There are B(M,K) possibilities to chose K elements from a group of M, where 
B (M, K) is the binomial coefficient. 

2. It is possible to have K! permutation of K elements. 

3. The escape probability of N elements from D detectors with efficiency e 

each, is given by: (l-De)**N. 

So the probability of getting K hits in an array of D detectors all with 
efficiency e, when we allow any order of detection, is given by. 


For calculating hit probabilities you need to know: 


1 . 


Wio i7*<- hi' I 

B(M,K) *K! *B(D,K) *e**K*SUM [one escapes, two escapes, three escapes etc. ] 

You will notice that the sum stands for single hit (i.e. M-K escape detection) 
double hit etc. etc. While we know that the probability of double hits in each 
one detector is very small, when summing over all the combinatorical 
combinations a large number arise. So when you do this algebra, as I FAX to 
you, you will recover our formula in the DOE report. 


Steve in our last phone conversation you promissed to me that once we solve the 
question of calculating the efficiency of our neutron detector, you will agree 
to publishing our results. Do you agree now? 

■far 


Best Regards Moshe Gai. 
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Response to Comments by Dr. Gai to My Analysis 
of the Yale Cold-Fusion Experiment 

A.N. Anderson PhD 
9 Jan, 1990 


Background 

This is a response to comments by Dr. M. Gai of Yale University, on a report I did as a 
consultant for Dr. S. E. Jones, a collaborator with Dr. Gai, in an experiment which, in 
view of the impact it may have on future research, he felt should be carefully scrutinized 
to see what conclusions I felt were justified. 

My report was to be only for the internal use of the authors, which they could use or 
not as they saw fit, and in that light I tried to be constructive. The questions I raised, I 
feel deserve real answers. 

I will address the comments in the order they appear in the communication from Dr. 

Gai: 

• Credentials: Dr. Gai assumes the analysis was done by a student of Dr. Jones, and 
as such is perhaps not to be taken seriously. Whether that opinion of my report is 
justified or not, my credentials are not irrelevant and I present them here. I received a 
PhD in Nuclear Physics from the University of Idaho in 1975. I have taken part in and 
analyzed experiments at LAMPF (Thesis on Proton-Proton Elastic Scattering, muon- 
catalyzed fusion), TRIUMF (University of Alberta, Proton-Proton Bremsstrahlung, 
neutron scattering, and muon-catalyzed fusion experiments) and Rutherford Appleton 
Lab, U.K., (muon-catalyzed fusion). 

Specifically I have been analyzing neutron-burst data from muon-catalyzed fusion 
experiments at LAMPF since 1983, in collaboration with Dr. Jones. The problems 
associated with these experiments are very similar to those encountered in the Yale 
experiment, and I can 6ay with 6ome assurance that I am as qualified as anyone to 
comment on them. 

• Was this an analysis? My report included, among other items, a plot of burst de- 
tections, and two tables of neutron singles counts. I was able to infer both neutron 
singles and burst rates (within efficiencies). What I was unable to do was develop 
a satisfactory feel for the actual efficiency of the system. I raised questions which I 
think should be answered before conclusions are drawn. 

• Flaws and TDC’s: It is not my function to discuss why flaws occurred in the exper- 
iment, but rather whether they occurred. Time of flight is an extremely important 
criterion by which background gammas, etc. can be rejected. The way the experiment 
was 6et up, however could not allow time of flight from more than one neutron per 
central detector, and in the requirement of the latch (of which more below) means that 
time of flight was unavailable for events occurring in the second center detector over 
500 ns. after an event in the first detector. I belabor this point because a little reflec- 
tion would have shown that high resolution was necessary only for gamma rejection 
over a very short interval, and was certainly not needed over the entire 20 /is. Dr. A. 
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J. Caffrey solved a similar problem for our experiment at LAMPF in 1983, using logic 
routers and groups of TDC’s of two different resolutions, which were readily available. 

• Crosstalk: It is irrelevant whether crosstalk is common (In fact, I have never observed 
TDC crosstalk at all, let alone this magnitude, even in the 50 TDC’s we have used for 
the last 7 years at LAMPF). The proper question is, does it cause problems? Figure 
2 in my report shows a small number of events in detector D5 (run 52, Background) 
which satisfy all criteria I could put on them except time of flight, (i.e. Both latches, 
Central=neutron, Ring=gamma). Besides those with good times of flight, there are 
several events whose times of flight look suspiciously like the gamma peaks from other 
detectors. If the peaks of different detectors had been lined up (as several others were) 
there would be no way to distinguish them from good events. It is the burden of the 
authors to demonstrate that these events do not compromise the data. 

• Latches: Dr. Gai implies that not only is it not necessary to look at latches, but that 
requiring latches throws away events! My plots of event gates show^ that- almost all 
otherwise good events also had latches. A very large percentage of events which did 
not have latches had problems (like leftover TDC’s and crosstalk). Requiring good 
ADC pulses does not allow the deadtime of the ADC to be evaluated. Likewise for 
the pulse shapes. The cause of bad pulses is likely to be pileup and as such reflects a 
decrease in efficiency. The effects I observed were as great as 70% in some detectors. 
The. authors have a responsibility to show whether this effect is real and if so, has 
been accounted for. 

• Saturation: As I stated in my report, amplifier saturation was not in itself a fatal 
problem. I was struck more by the lack of uniformity of all aspects of the pulses, 
including saturation, gain, thresholds and pedestals (Fig. 5, my report). Again, if I 
thought that these variations occurred in spite of the extreme care that had been 
exercised, I would not consider them important. It is the suspicion of sloppiness 
that worries me. When an experiment is acknowledged to be pushing the limits of 
experience, it is all the more necessary to take extreme care to ensure that these 
particular parameters, are carefully set and monitored. It is puzzling why some of 
the counters would have their pedestals carefully set to or below zero. Again, it is the 
burden of authorship to answer these questions to the satisfaction of readers. 

• Probabilities: This is a most exasperating point. When I first saw the expression used 
in the ‘Interim Report’, I assumed it was the work of a naive student. I am astounded 
to see the same absurd expression being defended by a professor of Physics at Yale. Let 
us quickly dispose of the arguments against my expression. The particular expression 
I used is indeed the probability for one and only one hit in a detector. There were 
two reasons I chose it. First, it was as applicable to the experiment as any other 
assumption about multiplicities. (The argument has been made that events with 
any multiple hits might need to be rejected. In fact, that was the stated reason for 
rejecting the single 5-fold hit mentioned in the ‘Interim Report’. See my report for 
discussion.) Second, it was easy to derive, and evaluate, and so was a good check on 
the MonteCarlo, which was my primary source of probabilities. 

I examined the source from which Dr. Gai claims he derived his expression, namely 
L. Westerberg et al -v-Rav Multiplicities From Nuclear Reaction s: (etc.) , NIM 145 
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(1977) 295-310. I was unable to find anything which remotely resembled the expression 
in the ‘Interim Report’. The closest I could find to a relevant expression was the 
following (with variables changed to try to correspond with my expression): 

Pm = (-1)‘ (f) £(-l)‘ (;) [1 - (D - •>-]" ni 1 - - •>/>) 

Let met point out explicitly, that this expression is intended to be the probability of 
one or more hits in k detectors. It in no way contradicts the expression in my report! 
The expression looks plausible to me, it satisfies the more obvious requirements of a 
probability, but having derived several such expressions over the years, I have learned 
not to trust such complicated expressions without testing. I am not sure about the 
relevance of this expression to the ‘Internal Report’ and will leave testing of it to. 
someone else. Let me dispose of Dr. Gai’s arguments instead. 

Consider Dr. Gsu’s examples. For 100 % efficiency and multiplicity 2 , the probability 
of detecting one and only'one neutron is indeed zero, and an expression which gave 
any other answer should be rejected. Again, for 50 % efficiency, one would expect 
my expression to give decreasing probabilities for large bursts, and again one would 
be right. This is exactly what one would observe if one rejected events with more 
than one hit. This is what will happen if one requires pulse information from central 
counters during large bursts. Let us now examine my expression for one detector 
of 50 % efficiency and a burst of 2 neutrons. Simple probability tells us that the 
probabilities for 0 , 1 and 2 hits are respectively, 0 . 25 , 0 . 5 , and 0 . 25 . My expression 
yeilds 0.25 and 0.5 for the first two, and is not defined for the third, since it would be 
the probability of having 2 detectors get one and only one hit. 

Dr. Gai claims that the 200 % his expression yields can be interpreted as a 100 % 
probability for a double hit. This is certainly not a probability or even an efficiency. 
It cant be an expectation for number of detectors hit, because it can easily exceed 
the number of detectors. Using this same reasoning, should the 100 % one gets from 
2 neutrons on a detector of 50 % be interpreted as a 50 % probability of a double hit? 
Perhaps a 33% probability of a triple hit? 

Since this is not a correspondence course in elementary probability theory, I will move 
on. This discussion should not, however, be dismissed, because the differences in the 
rates it predict are so significant. It definitely is not a trivial point! 

• Advertized Veto Efficiency: I use that word because, not having measured it myself, 
I was taking it as given. Perhaps the word ‘reported’ would be more appropriate. 

• Analysis hampered by lack of diagrams: As I was eventually able to infer many of 
the properties of the electronics from the data themselves, the only thing that was 
lost was time. The authors seem to think that time is a fairly important commodity, 
however. 

• Communication: I have spoken over the phone to Ralph Ranees and found him 
helpful, but he admitted to some ignorance of the electronics. I left several messages 
for Steve Rugari, but he did not return my calls. 
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» Detector Efficiency: I reiterate, that in any experiment that attempts to measure 
neutron rates, detector efficiency is of utmost importance, and should be completely 
understood before publication. I am somewhat puzzled by the rush, however. Is it 
more important that the ‘Interim Report’ be pubUshed soon than that the results be 

correct? 

Conclusion 

I suggest that my report be taken seriously. Most of the points mentioned by Dr. Gai could 
have been cleared up by careful reading of the report. Several of my worries and questions 
may already have very satisfactory answers, and others may not significantly change the 
results, but they are all legitimate. The question of the probability distributions should 
be seriously addressed. I suggest consulting a recognized expert on probability. If tlus 
problem cannot be corrected, I would advise Dr. Jones and anyone else connected with 
the experiment to remove his name from the author list to avoid the real possibility of 




To: Yale/BNL/BYU collaboration 

From: S. E. Jones 

Date: 11 January 1990 


Dear Colleagues: 

I find it an unpleasant task to respond to Moshe's comments on Alan 
Anderson ' s report, that "it is so bad that we might as well just 
ignore it," that "I don't see any plots of bursts detections, 
random emission etc. per run," that Alan "managed to confuse you 
on the issue of calculating hit probabilities," and so on. A1 has 
responded, and I agree with him that a careful reading of our 
reports would have cleared up many of the questions raised by 
Moshe . 


For example, random neutron emission is discussed on pages 1,2,3, 
9 (relevant plot) and 15 of Al's report in some detail, including 
a table (p. 2) showing the per-run neutron singles rates and an 
appendix (p. 15) showing the actual singles counts in each ring 
detector per run. Burst detections are discussed in my report, and 
in Al's on pages 3, 4, 5, 6, 7 and related plots /(Figures 3-7). I 
would like to call particular attention to Al's Fig. 4 which 


relates directly to the burst-candidate in 

locotflon of 3rd ring neut. 

‘ from Vole analysis 


run 


45: 



u#s 

45#C 

Figure 4. Pulse height spectrum lor U5 in 
data run 45, showing location of the triple- 
neutron coincidence event which was rejected 
for haring too laxge a pulse height. 
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Here A1 shows the location of the pulse height (PH) and pulse shape 
(PS) for the detected hit in detector in detector U5 (see also 
Table I in my report and Table III in the DOE interim report) . A1 
shows here that this point (PH=720, PS=1300) lies just above his 
2-dimensional "neutron gate" and also near the bottom edge of the 
saturation band for this detector. As I argue in my report, this 
point could represent a double hit in detector U5, resulting in an 
energy higher than 2.5 MeV (about 4 MeV). But the saturation band 
interferes with an unambiguous interpretation, since the point 
could be associated with the saturation band -- a gamma, for 
instance. The saturation may thus reduce our ability to discern 
double hits; this responds to a couple of Moshe 's questions. I 
suggest that the saturation band situation biases against high 
neutron multiplicities since these tend to have more multiple hits. 
The ability of the set-up generally to handle multiple hits is a 
problem I discuss also in my report. 


Moshe excuses the lack of relative timing between the two sets of 
detectors (up and down) as due to lack of sufficient TAC's. But 
since these detectors are much more efficient than the ring 
counters, it would seem logical to sacrifice a ring counter's 
timing in the event of too few electronics modules in order to get 
timing information from the two counter sets. This oversight 
greatly impedes our ability to extract useful timing information 
from neutron bursts, as I discuss in my report. It is a major 
weakness in the set-up, in my opinion. 


Moshe states that requiring pulse height information 
eliminates the cross-talk. However, A1 demonstrates in Figure 2 
of his report that evident cross-talk persists even with pulse- 
height information (and latches) required. (A1 discusses this 
further in his response to Moshe.) I find that A1 has done a good 
job of supporting his conclusions with tables and plots. _ A1 has 
had about seven years of experience in looking at data coming from 
detectors such as those used at Yale and as such is in an excellent 
position to help us understand the data. We should be grateful 
that he was willing to help. 

Now we come to what I argue in my report is a central issue: 
calculating hit probabilities. In his report, A1 clearly states 
that the hit probability equation used is just for the "simple" 
case of having K detectors hit by 1 neutron only: 


v0- 


. The proper expression, verified by 
MonteCarlo, for the probability pM,D,t(K) of haring K detectors hit by exactly 1 neutron, 


and no other hits, given a burst of M neutrons and D detectors of efficiency e, is 


Pm,d.< m = QQ - De) M - K 

binomial 


( 1 ) 


where (^) is the usual binomial coefficient. Other hit criteria or differing efficiencies 
can make the expression much more complicated. The report uses an expression which 
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To handle more complicated cases such as multiple hits (and 
detectors of differing efficiency, such as ring and central 
counters— see Al's report p. 3 bottom), A1 uses a Monte Carlo 
calculation, an excellent approach. He then checks the Monte Carlo 
with the expression for consistency in the case of single counter 
hits (compare the last column [predicted probability based on the 
formula above] of his Appendix B. tables with the third column 
[relevant Monte Carlo calculation]). The agreement is good 
(remember, too, that just 10,000 throws were used in the Monte 
Carlo calculation per case). 


Again, we rely on the Monte Carlo, not the formula, for 
multiple— hit cases. A1 clearly states this in his report. But 

Moshe says the formula is wrong because it dres not allow for 
mu 1 1 i p 1 e hits: 
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Notice that Moshe's exDression agrees with Al's for the case of 
single hits, which is, of course, all that Al's expression was used 
for. Moshe's expression for multiple hits, however, is 
questionable and there are definite errors in getting from that to 
the (incorrect) formula used in the Interim Report. A1 has 
explained these problems in detail and I will urge you to carefully 
read his comments. ^ o 
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Faced with the same ( non-trivial ! ) problems as described in 
our long reports, I draw the same conclusions on the best path to 
follow now in our quest to do a good experiment. I think we should 
take what we have learned and carefully design and execute an 
experiment to examine the time— structure of possible neutron burst 
events . I have recommended the use of the segmented neutron 
detector at LAMPF (Los Alamos) for the reasons given in my report. 
This detector has been successfully used in the past to analyze the 
time structure of bursts of neutrons arising from muon-catalyzed 
fusion and can be applied to the task at hand. 


However, I have learned to my dismay that someone called the 
director of LAMPF (G. Garvey) last Tuesday (Jan. 9, 1990) and urged 
him NOT to allow us to do this particular experiment! My response 
is that there are other places where we can do the work, if 
necessary, including foreign countries where we have collaborative 
arrangements with muon— catalyzed fusion research groups. 
Scientific research can clearly be impeded, but it cannot 
(ultimately) be blocked. 


Gentlemen, we intend to continue in this research and to do 
our best. I extend an invitation to all collaborators to come to 
BYU soon for a discussion of the problems and the best solutions. 

Sincerely, and perhaps a bit more 
firmly than usual for me, 

f 

Steven E. Joii 

cc: R. Gajewski, Douglas Morrison, Peter Parker (since Moshe has 

them in related communications , to which we are responding) 
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To: Steve Jones < j onesse@YVAX§ : 7^^Tv3H^^ynn <kgl%APSEDOFF .BITNET@VM>, Douglas 
Morrison <jeremy%UCBCM^ACBlTNET@VM> 

X- VMS -To: Steve Jones onesse@YVAX . BYU . EDUX kelvin lynn 

<kgl@APSEDOFF . BITNET^ , Douglas Morrison ^jbremy@UCBCMSA. BITNET> 

Resent-date: Thu, 11 J^n 90 15:32 
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Message-id : <F5EF8043355F2000BA@YVAX . BYU . EDU> 

X-Envelope- to : JONESSE@VANLAB . BYU . EDU 


Dear Steve, 




V' 


It was good to talk to you on the phone. I am happy that you disagree with the 
style of writing of your assistant A1 Anderson . 

The problem of calculating hit probability ijsT very s impler-arid exact solutions 
were given by B.R. Mottelson in 1974. In the interim report we wrote down 
first order approximation of the exact solution. It is O.K. for low 
multilicities as is the case here. 





From my graph you can see that for a multiplicity off 20/we have/ 85% chance of 
obsepying fold 4 or 5 or 6 or 7 or 8 or 9 or 10 . Sin^e- we do nqt^gbs er ved in 
the ([ring) counters more than a fold=3 event, than with 85% confidence we can say 
that no events of multiplicity larger than 20 were observed^/ h ? 


Best Regards Moshe Gai . 
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To: Ya 1 e/BNL/BYU collaborators and interested parties 

F r o to : S * E . Jones 

Date: 15 January 1 S 9 0 


Dear Colleagues, 

As you know, we have been carrying on a dialog regarding the 
interpretation of the data acquired using the Yale apparatus 
last August. Enclosed find the reports by myself and Dr. Alan 
Anderson, Moshe ' s responses, and our replies to Moshe. I would 
be willing to continue this dialog, except for two watershed 
factors which I learned about just last week and which I 
discuss below. The decision to discontinue working with Moshe 
has been made after careful deliberation this past weekend, and 
out of respect for all concerned, I provide the reasons here. 
Please understand that the attached documents provide relevant 
infer ma tion that I ask you to consider in addition to i a 
brief letter. 


I. The probability of detecting neutrons from a Menlove-type 

burst (approx. 20 to 300 neutrons) must be well-known in order 
to interpret the data, and indeed to assess the eff icacy of the 
experimental set-up. We have been using the following formula 

supplied by Moshe: lUW^dtfori'f^S) 

The most probable neutron Multiplicity (Mn) can be ccicuiecea 
from an eveht of fold (K) using the efficiency for detecting an event of 
fold K given by: 

€<K) - H n xe(12)*(M n -l)*«(ll)*...*CM n -K)* c(12-K) 

For neutron source of multiplicity 125 ve obtain 100% efficiency for 
detecting a neutron in each detector (fold-12), and for multiplicity of 20 
we have 100% efficiency of detecting 5 neutrons. Hence a 5 (4) foU event 
most probably corresponds to a source with neutron multiplicity of 20 

(24). .7". .11. 

Alan Anderson and I demonstrate that this formula is 
incorrect (attached reports). I concluded: "the design of the 

Yale experiment and the data interpretation in the Interim 
Report are to a large extent based on a faulty formula," and 
"the burst-detection efficiency given in the Interim Report is 
greatly over opt imist ic and the conclusions that follow from it 
are misleading." (Jones, 1/3/90 report) Here is Moshe's 

response (see more attached): _ — — - 
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Moshe, I felt that your stubborn defense of your formula, 

including justifying probabi 1 i t ies greater than unity, was 
obtuse, but still I thought we might reach an understanding. But 
then late last week (Friday) I received a copy of a letter to you 

from Richard Gar-win, dated November 1 . .1989. Dr. Garwin states 

that the neutron detection efficiency formula in the Interim 

Reoort is incorrect, and aclcs: 



"I would like to receive from you a corrected version of the 
manuscript which does not have this misleading statement, and 
which does the analysis right, 11 In your response, you agreed. 
These facts (which I learned late last week) lead me to raise 
strong objections. 

Why was I not informed of Dr, Garwin’s insight earlier? When 
Alan and I independently dug out the same conclusion, why your 
stubborn struggle to maintain the validity of the repudiated 
expression (see attachments)? Why argue that Alan Anderson 
"managed to confuse you [Jones] on the issue of calculating hit 
probabilities" in view of the insight from Dr. Garwin? 
Finally, how could you let the paper reside with the DoE/ERAS 
committee uncorrected, when you had learned from Garwin that 
the Interim Report was misleading? I note that your infamous 
press release which accompanied the release of the Interim 
Report to the media was at this same time that Garwin caught 
the error -- you should not then have gone to the press! 


In view of these matters, I am not willing to continue the 
collaboration with Dr. Ga i . But there is more. 

2. Mo she will not concede the seriousness of electronics and 
other problems in the Yale set-up which are plainly visible in 
the data analysis (see esp. Anderson report). Anderson and I 
show that the saturation in pulse-height information in some 
counters and the cross-talk between timing channels (which does 
persist when ADC information is required, incidentally) and 
other electronics untidiness cause difficulties in the data 
analysis. The lack of timing between up and down detector sets 
is a major weakness of the set-up, in my view (attachments). 

In my January 3 write-up, I noted these problems along 
with the low efficiency for burst detection of the Yale 
set-up. 1 proposed that we take what we had learned and do an 
improved follow-up experiment, either at Yale or at LAMPF . I 
then discussed the improvements which could be realized with 
the LAMPF neutron detector system, used in the past for 
muon-catalyzed fusion studies. I hoped that this approach 
would avoid embarrassment -- it is not unusual tc have to re-do 
an experiment to correct problems. 

But instead of accepting this solution or even discussing 
it with me, Dr. Ga i called the director of LAMPF (G. Garvey) 
last week I am not privy to the details of that conversation, 
but I was informed that "pressures" could stop the experiment, 
and that the possible reaction of the physics community would 
have to be evaluated before the experiment could be pursued, as 
a direct result of the conversation. Mind you, the proposed 
experiment would use equipment which is dormant while the LAMPF 
beam is off (until May) and had received prior informal 
approval from Jim Bradbury. It would be conducted with very 
few demands on LAMPF personnel. Certainly, Moshe’s phone call 
did nothing to help the follow-up experiment to which ha had 
been invited, and is not indicative of a cooperative 
col laborat or . 


I have explained salient reasons for my decision to 
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discontinue collaborating with Hoshe. Please accept my 
statement that I bear no personal animosity towards you, Moshs, 
although several of your remarks and actions in this 
collaboration have cut me deeply as I have told you before. I 
am saddened by these events, but hopeful that we can make use 
of what we have learned in the effort to better the research. 
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Richard L. Gar win 
IBM Research Division 
Thomas J. Watson Research Center 
P.O. Box 218 

Yorktown Heights, NY 10598 
(914) 945-2555 

FAX: (914) 945-4419, Telex: 137456 IBMRESRCH YKHG UD, 

BITNET: RLG2 at IBM.COM 


November 1, 1989 
(Via BITNET to GAI at YALEVM) 


Professor Moshe Gai 
A.W. Wright Nuclear 
Structure Laboratory 
Yale University 
New Haven, CT 06511 


Dear Moshe, 


,< 3 > 

I have seen your "Upper Limits on Neutron Bursts ..." but 
before I try to follow the analysis in detail, I wanted to 
ask you about something that caught ray eye. You say "For 
neutron source of multiplicity 125 we obtain 100% efficiency 
for detecting a neutron in each detector (fold-12), and for 
multiplicity of 20 we have 100% efficiency of detecting 5 
neutrons. 1 * I understand that with a 0.8% detector 
efficiency and the multiplicity of 125, you have an expected 
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number of neutrons counted in each counter of 1.000, but 
that by no means corresponds to "100% efficiency." In fact, 
it corresponds to (1-1/e) or 63% efficiency. 

Of course you know this, but I would like to receive from 
you a corrected version of the manuscript which does not 
have this misleading statement, and which does the analysis 

rl9ht - 

Thank you for providing this beautiful experiment. I hope 
that the further projected analysis is available soon. 

Sincerely yours. 


Richard L. Garwin 
Forwarded in his absence 
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RLG: jah:305%MG: 110189.. MG 

} RLG2 YKTVMV 11/01/89 . . , 

~ R.L. Garwin GAI8YALEVM 11/01/89 Interpretation of expt wath J 


Date: 

From: 

Subject : 

To: 

In-Reply-To : 


Wed, 01 Nov 89 16:51:57 EST 
Moshe Gai <GAI(?YALEVM> 

Re: Interpretation of expt with Jones, et al . 

RLG26WATSON 

Message of 1 November 1989, 16:19:20 EST from <RLG2@YKTVKV> 


Deaf Dr. Garwin, 


Yes I agree. Please note that the manuscript was written in a day. On October 
3rd I suggested that we submitt a report but Jones refused to submit a report 
to DOE/ERAB panel. I did not plan to submit one due to his objection. We only 
received Jones OK for submitting a report on the data which was analyzed over a 
month ago, by October 24th. and a report had to be written, Figures made, all 
reach an agreement and report be submitted by Oct. 27th 1 While the report 
includes data that we all agree on without mistakes it may have a few omissions 
as the one you mentioned. We plan to work on these deteails when Jones returns 
from Japan next Monday and submit a fianl report to the comitt ee . 

Thanks for your comment and with best regards Moshe Gai. 

- Moshe Gai RLG2@ WATSON ll/01/89*Interpretation of expt with J 
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In the example we have above with one detector with efficiency of 100% and a 
neutron source of multilicity 2, using my equation yield an efficiency of 200%. 
In this sense it is sinply interpreted as 100% probability for a double hit. 

In fact you will notice that the formula we wrote in our report corresponds to 
t “ e first binomial times K!*e**K, and the rest of the factors in (1) need to be 
rewritten and summed over . While having efficiency larger than 100% as 
suggested in my formula are indeed strange, they simply have to be 
explained (i . e . they are not wrong) . 

May I suggest that we discuss it so as to avoid all this waist of time on 
trivia . 



Notes on the ^ 

Derivation of the probability expression in the Yale ‘Interim Report’ 

A. N. Anderson 
11 Jan 1990 


Introduction , , . . ..... 

The derivation of the expression used in the Yale ‘Interim Report , or ® P*° “ 1 

P MD (k) of getting k hits among D detectors of efficiency e, when exposed to a bu 

M neutrons, begins from the following expression: 




( 1 ) 


where e‘ is the probability that k of the M neutron 

the probabilities that the remaining M -k neutrons hit exactly P V 

detectors, and none of the others. The expression used for A is 


A-(i -D‘T 


M-k-i 


( 2 ) 


by which the ith term is meant to be the probability that M - k - * of the remaining 

neutrons miss all D detectors. _ , , . 

The summation is then evaluated as a geometric senes yielding 


PM,D( k ) 


-too 


M-k 


k\e 


1-(1 -De) 

1 - (1 - De) 


(3) 


The approximation 
is invoked to yield 


(4) 

(5) 


(1 — x) M « 1 — Mx 

which is essentially the expression used in the report. Actually ‘^XTh^mtathe 
should be (M - k + l)/fc! but that is a minor point, considering all the other 


derivation. 
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Discussion 

The starting point of the derivation is correct provided the A’s are chosen properly. Re- 
quiring one and only one hit per detector replaces the summation with just its zeroth term, 
Ao = (1 — De) M ~ k , which is the probability that all remaining M — k neutrons miss all D 
detectors. This results in 




which is the expression I used to check the MonteCarlo output in my report. It satisfies 
the requirements of probability in that it is always between 0 and 1 and it sums over k to 
a number between 0 and 1, namely the probability for having an event with no multipie 
hits. 

It is easy to 6how, however, that equation (3) is not a valid probability. Simply 
let M become very large. There are only two terms containing M. namely and 
1— (1— De) M ~ k+1 . The first term increases roughly as and the second term approaches 
1, since 1 — De is less than 1 and its power goes to zero. Therefore the expression can 
become greater than one and is not a probability. 

One reason (3) is not a probability is that the terms do not contain factors for all the 
neutrons. To allow for other hits on detectors, the expression adds conditional probabilities 
that M — k — i of the remaining neutrons also miss all D detectors. Unfortunately, that 
says nothing about whether the remaining i neutrons hit a detector. If they fail to hit a 
detector, they have already been included in the previous term and we are counting them 

again. X\ v _ 

In order for probabilities to be addable, they must be disjoint. The ith term, however, 
is included in the i — 1st term if the ith neutron does not hit a detector. In order to make 
them disjoint, then, the t neutrons must hit a detector, and the ith term should be of the 
form 

A { cx e'(l - De)**-*"*' 


;tors for eau 


where now there are factors for each neutron. 

The magnitude of the error can be appreciated by the size of the correction, that is, 
each term has now been decreased by a factor of e', which when summed over a large 
M, can make an enormous difference. It also means that the sum in (2) is no longer a 
geometric series, and hence the form of (3) is no longer implied. The new form has a much 
better chance of satisfying the probability conditions. 

We are not done yet, however, because we have not counted up all the ways the i 
neutrons can hit the detectors, that is, we have not determined the coefficients of A,-. In 
general, there is much subtlety involved and it is very easy to make errors in determining 
these, and I do not choose to pursue this path further. Many seemingly plausible arguments 
have generated expressions which, when evaluated, normalized to more than 1, or which, 
when compared to a MonteCarlo for the same system, gave discrepant results. There is no 
reason to believe that the most simple system might not have an extremely complicated 
expression for the exact probability. It is at this juncture that a MonteCarlo, properly 
tested against simple cases, justifies greater confidence than a derived expression. 
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For the sake of completeness and to underscore my remark about subtlety, let me 
introduce an expression which has the advantages of legality, plausibility, and agreement 
with the MonteCarlo. It bears little superficial resemlancc to (1), however: 


p M , D (k) = Q(1 - (1 - - (D - m M - 1 (•> 


The second term represents the probability that each of k detectors was not mtssed by aU 
M-k + 1 neutrons. (The other k - 1 neutrons presumably having been required to hit the 
other k - 1 detectors. The last term is the probability that the M-k neutrons not required 
for hits were precluded from hitting the remaining detectors. It has the nice property that 
as M becomes large, all probabilities go to zero except for k = D whrdi yields 1, as it 
should. Presumably a proper choice of A’s would yield this result. It is difficult to prove 
that this expression is exactly correct, however, or that the correct expression might no 


be much more complicated. , 

Having shown that (3) is not a correct expression, I will mention a-so that approxi- 
mation (4) is incorrect. The easiest way to see this is again to make Af large. Whereas 
(1 - De) is always between 0 and 1, and raising to the power M-K + 1 preserves this 
property, note that the approximation 1 — (M — k + l)De is unbounded, and expression 
(5) becomes greater than one even for relatively small M. 

It is necessary to dispose of just one more erroneous proposition, namely, that somehow 
inclusion of the (M - k)/k\ term restores legitimacy to (5). As I have already pointed out 
(5) is invalid even including the term in question. An example has been used m which 
D- 1 and e is 100%. In this case k is constrained to be either 0 or 1 (It must be less than 
or equal to D), so Jfc! = 1. In this case the term takes on the values 2 and 1 for M - 2, 
and 3 and 2 for Af = 3. The term actually increases already absurdly large probabilities . 
It should be obvious that this term does nothing to help (5). 


Conclusion 

It is not unreasonable that anyone might err in trying to derive an expression for the 
number of unconstrained detector hits. It is a somewhat subtle problem. What w unrea- 
sonable is that a Professor of Physics at Yale, who is ostensibly expert at experiments m 
which problems of this type occur regularly, should have so little understanding o e emen- 
tary probability that not only does he fail to apply simple tests to his results to ensure 
their validity (or at least their plausibility), but that he does not even recognize such tests 
when someone else applies them, and even goes so far as to invent new interpretations of 


probability. 

I hope enough has been demonstrated that there 
point further. 


will be no necessity to pursue this 
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Dear Steve, 


* 


Vi 


\€ 


**r 


Thanks for the memo. Good that we agree that the time for Yale-BYU-BNL 
experiment should be increased from 77 hours to 103. We are making progress, 
even if it is slow. I would still request the following calrified: 

1. Since you talk about cylinder *hour, should you multiply the Yale time by 

4 to account for the 4 cylinders? _ * o £ 

/ * 

2. If indeed at Los Alamos you spent a long time for cracking the surface of 
the oxyde layer, should you not multiply the time we spent on the treated T1662 
by the appropriate factor to allow for the fact that these samples in our 
experiment are already pretreated and the deuterided already formed? 




H0 




h ry y- £ 2 O 


3. In the Menlove paper in Table III, I count 27 bursts, and you told the DOE 
that you have a grand total of 37 bursts. So in /fact the data in the Menlove 
paper is about 75% of your data. We than can discuss these published materia] 

For the Menlove data I find: (J i 

/'*y' i ^ ^ 

A. All cylinder bursts during cycle/5 lor smaller (i.e. before^ three days of 
running) . Where is the evidence thau-one needs to wait for a TungVtime? You 
even have 11 bursts at cycles 1, 2, 3 and 4, i.e. within a day or two. The 
Ti662 sample bupsts within cycle 3, 4! ! Again where is the evidence that we 
need to wait w^eks?? 


Z *4 ± 2 <? 
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B. S irjCjQ-^fche-J^en 1 ove resull 
use&. 4 cylinders, twice, i.el 
obserfevd., — byjj^V _ j~f 
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■s ugg e st 30-50% chance per cylinder and since we 7 JV'/ 


bursts/ should have been 
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C. If you insist that no background is included in the 13,000 hours than k “4 

before I agree to that number I would require a detailed list of all experiment 
performed. Before I can calculate this number by myself, I will not put it in 
our paper. During the first period of the Menlove experiment you rotated 
between a dummy and a real cell, i.e. measured background half of the time. So 
you either should not multiply this time by 4 (four cylinders) or appropriately 
subtract the background time . 

\yn - ()» 3 

I just made the following observation: 

While the Menlove paper covers a period of two months, it includes ''75% of the 



bursts-. Than for the rest of the five months you only had 10 bursts, i.e. 2. 
of the grand total? Steve it does not seem to me that you have improved the 
situation, in fact I would say that it got worst... / ,, j ^ 


Best Regards Mo she Gai . 



Jar.. 4, 1989 


Dear Mo she, 

I have returned from a family trip for the holidays and can 
address your questions. I should also say that our report on the 
analysis will be sent tomorrow by express carrier. It raises a 
number of questions which I hope you will address. 

I will answer your questions point by point. 

1. No, the normalization is done by mass, not number of 
cylinders. That is, we ran at Yale for 103 hours with about £00 □ 
of material per hour. At Los Alamos, the 13,000 hours were 
estimated based on the number of cylinders, each containing about 
100 g of material. So the appropriate scaling is the ratio of 
<103 X £00 ) / <13,000 X 100), as I have said before. 

£# None of the Ti662 samples used at Yale were “pretreated and 
the deuteride already formed. “ When we loaded the deuterium gas, 
the pressure dropped very little showing that very little 
deuteriding was occurring. <The used cathode materials did 
rapidly form a deuteride at Yale, but not the TiS6£ samples.) 

3. Your count of £7 bursts in the original Men love paper 
includes 6 bursts from electrolysis experiments as well as bursts 
smaller than £0 neutrons. My count of 39 bursts is for gas- 
loading experiments only and for bursts of £0 or more neutrons. 
The Men love Update shows the sizes of bursts in Table II, so that 
in the data set relevant to the submission to Nature, I count 1£ 
bursts of £0 or more neutrons. So in fact the Men love original 
paper represents about 30% <12/39) of the data which is 
consistent with the time periods involved <this answers your 
concluding remark. ) 

3A. You state that “For the Mer.Iove data I find: A. All 
cylinder bursts during cycle 5 or smaller <i.e. before three days 
of running)." However, the original Mer.Iove paper in fact shows 
bursts through cycle 9 and a typical cycle at LANL requires a 
full day so this would correspond to roughly 9 days of running 
(where did you get “before three days"?). Some runs do produce 
bursts during the first cycle or two. But stat 1st ical ly, we find 
as mentioned in the Mer.Iove October updates “A time interval of 
several days is [often! required between the gas loading and the 
first neutron yield from the LN cycles." Thus, at LANL running 
for 14 cycles is not uncommon. 

3B. About 30-50'/« chance per cylinder depends on long measurement 
times as I have said before. Also, to see 2— 4 bursts at Yale 
assumes burst — det ect i on efficiency comparable to that of LANL. 

0 u r analysis s h o w s t h a t t h i s :i. s n o t a t a 1 1 t h e ca sf t h is is 1 - 

the reoort in detail. This low efficiency is the rnaior or oh 1 e" 

1 n the Yale set — u p i n m y j u d g m e n t , rat h e r t h a n t h e s h o r t r u n n i n z 
time, as we must discuss when you study our report carefully. 



3C„ We could ask Ed Garcia for a more detailed list of "all 
experiment performed" as you now say* I think the estimate of 
13,000 hours suffices when the revised efficiency of the Yale 
set-up is taken into consideration — let’s discuss this when you 
have read the report- We scale by the sample mass rather than by 
the number of cylinders as discussed above, and the background 
has been handled appropriately (i«e., not included in the count 
of 13,000 foreground hours) « 


Si ncerel y, 
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\v 


Dear Steve, 

O' . "x 

Thanks for your FAX note of Jan. 4th. Before replying to iir'let mb inform you 
that I talked to Dave Lindley and he is pretty upset that /you claimj that he 
misquoted you in the Nature note on Nov. 9th. I would then like to ! have a 
conference call involving all three of us so that we can get to the bottom of 


this . 


concerning your note: 




> 


jr r/>v / 


point 1: While for (DOE and to me(you quoted a burst chance per cylinder, you 

now normalize by weight? Which method are you eventually going to use? In the 
Menlove paper I see an average weight of 55 grams per cylinder, why did you use 
100 grams? For the Yale data you use real time 103 hours, but for the Los 
Alamos data you multiply the hpups-by the number of cylinders used in the four 
detector systems. We had our( four cylinders in one detector system, and again 
why have you omitted the factor of 4? I am afraid I will have to insist on 
receiving a detailed account of~“aT]j experiments performed in Los Alamos 
including: cylinder weights, time, /etc. before I agree to the 13,000 
cylinder *hours estimate. ‘ — ft () ifepett) 

point 2: I shall remind you that on October 31st I FAXed to you comments on 
our paper by Kurt Zilm. I felt that the comments directly criticized your 
experiment at Los Alamos and thus should have your approval . Since we did not 
hear from you they are still not included in the improvements to the report. 

The comment read: "Before the Ti662 alloys used in the majority of these 
experiments were loaded into the pressure cells they were cleaned according to 
the protocol developed at Los Alamos . This was done to ensure that our samples 
were as close as possible in composition to those investigated in Ref. 2 
(Menlove) . Since the oxide layer on the alloys is not removed in this 
treatment it is higly unlikely that a significant amount of deuterium is 
incorporate into the metal lattice. . and as you write in your point 2, Kurt 
Zilm adds: "The sample in set E which contained used Pd electrodes did have a 
significant amount of deuterium...". So in fact what was originally a severe 1 
criticism of your experiment at Los Alamos (by our material science people) is / 
now turn around by you as a cricism of the Yale experiment? . ' 

fh<i ijojf of -- 5 ca u ’ u c 6 r ret. fad 

point 3: It is true that 6 bursts are from electrolysis samples. But I shall 

repeat, in your paper with Menlove you explicitly write: "the neutron emission 

always ceased after a few cycles"; and before that: " A given cylinder was put 
through 7 to 14 of these cycles" . Each cycle lasted about a day with some 5 

a bit smaller than the number of cycles) . We 


hours (so the number of days is 
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hot di S a yfjVf ^ ) 

in fact run 12 cycles, just like Menlove. Do you still stand by the above 
statement of your paper? U/ to ^ zHo^s 

point 3B: Let me state for the record you are now claiming that: "Also, to see 

2-4 bursts at Yale assumes burst-detection efficiency comparable to that of 

LANL. Our analysis shows that this is not at all the case-- this in the report 

in detail. This low efficiency is the major problem in the Yale experiment 

set-up in my judgement, rather that the short running time, as we must discuss 

when you study our report carefully". Let us make it very clear you are now \) fe -f. 

taking a new path in your attempt to kill the Yale experiment. I am looking 

forward to receiving your study of the burst detection efficiency of the Yale J ' , 

exper iment . I hope you realize that calculating the efficiency is detailed, 

and just make sure you include all the factors before you mail to me your . j 

report. In fact you should consult our report where we quote the formula for J'^~ f 

the efficiency on p. 5. I noticed that Menlove in his letter to me made a 
mistake in calculating the efficiency of our experiment, and just make sure 

that you don't make the same mistake. . . /„ 

1 ooch> 

With Best Regards Moshe Gai . /, i • 




From: IN%"RLG2%WATS0N . BITNET@VM" 19-JAN-1990 09:52:25.68 

To r JONESSE@VANLAB . BYU . EDU 

CC: 

Subj : Retransmission. 

Received: from YVAX.BYU.EDU by VANLAB.BYU.EDU; Fri, 19 Jan 90 07:53 MST 
Received: from JNET -DAEMON by YVAX.BYU.EDU; Fri, 19 Jan 90 07:52 MST 
Received: From YKTVMV (RLG2) by BYUVAX with Jnet id 8416 for JONESSE@BYUVAX; 

Fri, 19 Jan 90 07:52 MST 
Resent-date: Fri, 19 Jan 90 07:52 MST 
Date: 19 January 1990, 06:59:03 EST 
From: RLG2%YKTVMV . BITNET@VM 
Subject: Retransmission. 

Resent-to: JONESSE@VANLAB.BYU.EDU 

To: JONESSE@YVAX 

Reply -to : RLG2%WATSON . BITNET@VM 

Resent -message- id : <EFE66D178F3F200733@YVAX . BYU . EDU> 

Message-id: <EFE678E9609F2000BA@YVAX.BYU.EDU> 

X-Envelope-to : JONESSE 

X- VMS -To : JONESSE@YVAX . BYU . EDU 


Date: 

From: 


To: 

cc: 


18 January 1990, 22:07:07 EST 
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Subject: NOTEs of 01/18/90. 
Reply-To: RLG2 at WATSON 


& 


Thanks for the NOTES, received fine via RLG2 at WATSON. I was most interested 
in your letter to Howard Menlove, and dismayed that you have not yet received 
a response. 

Certainly his criticisms of (a) your "low efficiency" and (b) your inadequate ~)~) 
running time are unwarranted. I will be most interested in any data Howard 
has obtained (ALL data, in fact) while running in caves or under other 
circumstances . 


V. 


Indeed, when I pointed out that your formula could not be "CORRECT" 

I did not mean that it was really misleading in the range where it would be 
used-- that is, for efficiencies substantially less than 50% or so. And it 
is never VERY far off. 

I would be grateful if you would put copies of the references into the mail 
for me-- it will save me a good deal of time in the library. 

Thanks VERY much for the information. 


V 1 

H*** 1 


ds 


Dick Garwin 
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Please note that Jones' e-mail is JONESSE @ BYUVAX (your cc has JONESS...) 
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"MOSHE GAI , (203)432 5195, FAX: (203)432 3522" 


From: IN%"GAI%YALEVM . BITNET@VM" 

To : J ONES SE@VANLAB . BYU . EDU 

CC: 

Subj : 


Dear Kelvin and Steve, 

I have written draft number 2 of our paper that includes answers to all BYU's }{cui ! j? 
questions, including the exact solution of hit probabilities and the various 
approximation that can be used. Note that in draft two I only take the exact 
Mottelson solution, without any multiple scattering (one neutron hits two 
detectors) or double hits (two neutrons hit one detector) . This was done in 
draft number 1. Since I feel we can easily agree on Mottelson 's formalism, and 
in order to expedite publication I am willing to stick to the most simple case, 
even though multiplee scattering and double hits make our efficiencies even 
larger. Remember the probability for such events in each detector are small, 
but when you work out the combinatoric for the whole array, it makes a 
difference. I have FAXed a copy of these calculations to BYU last week and I 
am EAXing one to you Kelvin now. 

I would be happy to add any more answers to unresolved questions in future 
drafts. I would think that the differences between BYU and Yale are now so 
small that the whole discussion can be resolved within the collaboration, and 
we should proceed to publication. I will hgweveis^Tiave no objection to 
arbitration! I would only hope that the^^arrasmeriE^paiised by the letter 
from Dick Garwin (cc to Jones) would bring tne voice - of reason. 

Best Regards Moshe gai . ( j/ni/ill'tf n nJtr/^ J 
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From: 
To: 
CC: 
Subj : 


Cleaner analysis of the data 


Dear Steve and Kelvin, 


BRIGHAM YOUNG 
UNIVERSITY 

THE GLORY OF GOD 


In the data analysis for the interim used only one dimensional gates 

on the pulse shape spectra. Ralph has also performed the analysis with gates 
on the two dimensional pulse shape spectrh. *1 discussed these results with you 
Steve over a month ago, and ('we both know thatythis more conservative analysis 
yields an amazingly clean spWtrjaw Jfe-are- how going to incorporate in Draft 
numbers— 2--olL_our paper the results c|f this analysis, too. The results are 


embarrasingly cleanN even cleaner 
Bestr-Reqard s M gshe Gai . 


an what we have gotten so far ! ! 
.UstJoi t jji)" t/bovt ft j G 
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January 18, 1990 

Dear Moshe, 


Quoting authorities (G. Hagemann, 0. Andersen, L. Westerberg) to 
justify your hit probability estimates is not sufficient -- your 
work must be correct. You mention that your formula is an 
approximation which immediately should raise worries, but A1 
shows in detail where. you went wrong <11 Jan. memo) -- please 
study and respond. I am truly amazed to read in your memorandum 
of yesterday: "I shall state again, there is nothing wrong with 
Moshe's formula as appeared in the [interim] report." 

In the latest Report, you state "we estimate for multiplicity M > 
22, at least 92% hit probability of 4 or more neutrons in the 
ring counters (i.e., P4(22)+... +pl0(22) > 92%). Hence, from our 
data we deduce an upper limit on neutron bursts from our cells, 
to include less than 22 neutrons with 90% confidence. " The 
conclusion is terribly misleading since tha_ hit probability on 
which it depends is way wrong (just as in the original Interim 
Report) . 

I wish you would respond to Al's detailed discussions and Monte 
Carlo calculations regarding hit probabilities. I have 
summarized in Fig. 1 the hit probabilities for 4 or more detected 
neutrons in the ring counters [ P4 <M) ] . Notice, for example, 
that the probability P4 for M = 23 neutrons is about 5%, not 92% as 
you obtained. That is an enormous difference and drastically 
affects conclusions obtainable from the experiment. Al's 
probability for one or more ring counters to be hit is about 80% 
(Al's tables for M=20 actually), which is still less than your 
92% for four hits. 

Is your value reasonable? A simple calculation suffices to show 
that it is not: 

23 source neutrons X 0.0008 (efficiency/ring counter) 

X 10 ring counters 

= 1,8 expected hits. 

But then the probability of detecting 4 or more neutrons must be 
considerably less than the 92% quoted in your latest report! 
Your approximate formula evidently breaks down even for small 
multiplicities like this and thus has lead to misleading 
conclusions as we argued at length in our Jan- 3 and subsequent 
reports. (A1 even shows [Jan. 11 memo] where your expression 
double-counts neutrons; you have not responded.) You will find 
that the expectation value of about 1.6 hits for 20 source 
neutrons comes out of Al's Monte Carlo calculation just as 
expected (Al's report of 6 Jan., Appendix B) . 

The fact that your formula checks for M=1 (as quoted in the 
latest report) is hardly convincing -- the probabilities are 
already far too high for M=23! 


In addition to this gross error 


the quoted efficiencies of 0.8% 


per ring counter are suspect. We have run a Monte Carlo for the 
approximate geometry of the Yale set-up and get signif icantly 
smaller values. (Remember, if you measure the efficiency with a 
source, a simple one-dimensional pulse-height/pulse-shape gate on 
neutrons will lead to inaccurate results.) A1 shows in his Jan. 
6 report that serious electronic i nef f iciencies are evident 
(e.g., counters U3 and D1 ) . For these reasons, I have included 
on Fig. 1 the probabilities obtained assuming ring detector 
efficiencies of 0.5% which may still be too high. Notice how 
insensitive the set-up is to bursts of even 50 neutrons! 

I find that my challenge of the hit probability statements in the 
Interim Report and the latest version corresponds quite closely 
with Richard Garwin's challenge (Nov. 1 letter to you, which I 
learned about only a few days ago) . And I will ask just as he 
did then! 

"I would like to receive from you a corrected version of the 
manuscript which does not have this misleading statement, and 
which does the analysis right." 

Be sure to fold in the correct geometrical and electronic 
efficiencies for each counter. And include in the analysis the 
(poor) handling of multiple hits in counters in the Yale set-up- 
- the Monte Carlo calculations show that these are not of low 
probability as erroneously stated in both the Interim and latest 
reports. Please do not ignore the observed excess of hits in two 
ring counters in foreground data versus background (Al's Table 1, 
Jan. 6 report). Finally, when comparing with Menlove, be sure to 
include the fact that the Yale set-up is sensitive to bursts for 
only 20 jus, while Menlove is sensitive for up to 128 ^s. Bursts 
could be quite long and spread out — a recent burst candidate of 
170 neutrons detected here lasted about 150 ji s. 

A1 and I have provided lengthy documents which attempt to do the 
analysis right, and we feel that your responses do not adequately 
address the concerns we have legitimately raised. This is 
particularly true of our hit probability calculations . The 
latest draft of the Interim Report incorporates nothing of our 
analyses in it that I can see. This is not acceptable. 

I find myself repeating, repeating, repeating what I have said in 
past memos and reports. Look, I really think the best thing to 
do is to admit that the Yale set-up was far from optimal for 
neutron burst detection and to do the experiment right based on 
we have learned, THEN publish. We can make an apparatus MUCH 
more sensitive for burst detection and do the definitive job we 
set out to do (whether we see bursts or not) . In my January 3 
letter to the collaboration, I suggested this solution. I think 
it is far better than rushing to publish based on the set-up of 
last August. 


With best wishes to a stubborn associate. 
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Pl(l) -the hit probability, for one neutron in any'~ofig- uf u he\ ten rine 
counters, is 8% , as ''expectec^, We estimate for multiplicity M>22, aV 1 
92% hit probability of 4 or more neutrons in the fringjcounters (i.e P^(22)+ 
... +PloC22) > 92%). Hence, from our data we deduce 


an upper limit on 

/ ^ 

neutron' bursts from our cells, to include less than 22 neutrons wicti 90% 
confidence. Using the data of Menlove et al. [2], we expect some 2-6 

neutron Wsts, including more than 20 neutrons each, over the duration of 
our experiment. - 

In Fig. 6 we show the rate above background of neutrons emitted 
randomly. For these data the time of flight coincidence method was used 
[3], that yielded a low background of 2 cph (0.4 C ph in upper hemisphere 


Laftit \y irs ro - of ^<50o/ft 

ni QS kl?) 

event. Since the probability for a double hit in a ring counter is 

estimated to be small er^than 3%, ye conclude that all events are consistent 
with background events, with 97% confidence. 

Inspecting the time scale and energy deposited in detectors U5, U1 
and U0 in the event of Run 45, we may suggest that the event is consistent 
with one high energy (E > 7.3 MeV) cosmic-neutron (arriving from the less 
shielded front side, see Fig. 1), that undergoes triple scattering, 
including one large angle scattering from detector U5 to its immediate 
neighbor detector Ul, and then followed by a small angle scattering into 

the central detector U0 lying below the ring detectors (thus yielding to 
negative time of flight) . 

We have calculated the hit probability in our array using the 
formalism first developed by B.R. Mottelson and discussed in Ref(s) 4 The 

u At 

hit probability is essentially determined % the array of ring detectors 
In Fig. 5 we show the hit probability -P K (M) , for the ring counters only, 
as a function of multiplicity (M) , for different-foids^IO - Note that 
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COLLEGE OF ARTS AND SCIENCE 

DEPARTMENT OF CHEMISTRY 


IA 


January 3, 1990 


Dr. Richard L. Garwin 
IBM Research Division 
Thomas J . Watson Research Center 
P.0. Box 218 

Yorktown Heights, NY 10598 
Dear Dick: 


A 


0 




I appreciate that you have sent me several relevant pieces of information 
on cold fusion since the last meeting of our Panel on October 31. Since I'm 
attempting to keep an up-to-date file on cold fusion, I would appreciate any 
new information that comes across your desk. 

rfO 

The final report of the Panel has now been duplicated and you should get a 
copy soon. It took some three weeks beyond October 31 to remove all the 
typographical errors and make minor changes, however, most of the delay has 
been due to DOE's reproduction schedule. 

Happy New Year. 


JRH : dsm 
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TEXAS A&M UNIVERSITY 

DEPARTMENT OF CHEMISTRY 

COLLEGE STATION, TEXAS 77843-3255 

December 26, 1989 


(409) 845-5335 
FAX(409) 845-4205 


Dr. John Huizenga 
Co-Chairman 

Professor of Nuclear Engineering 
Department of Chemistry 
University of Rochester 
Rochester, NY 14627 

Dear Dr. Huizenga: 




0 


I am now writing to you in respect to the present position of 

experimental work on cold fusion and its connection to the JrRAB report on that 
subject. 


With the recent publications describing/nucrfear electrochemical effects 
from Oak Ridge, Los Alamos, Brookhaven, and(N.R.C/‘ the Fleischmann-Pons 
phenomena have now been repeatedly conf irmed<by-tf. S . Government Research 
workers. These men have added their confirmatory experiments to those 
obtained from April onwards at Texas A&M, Stanford, Minnesota and elsewhere. 

part from the U.S. Government labs and universities, these startling 
phenomena have been observed and reported on in many other countries. 

As to those expert and respected scientists (e.g.. , at Yale, MIT Cal 
Tech, etc.) who denied the existence of these effects, it is now easy to 
understand the absence of their results: the effects, now seen by so many do 
not switch on for many weeks of electrolysis. As to reproducibility it is 
poor. However, for 1 mm wires, we have a 70% chance of producing tritium if 
electrolysis is continued for more than 10 weeks.. 

as now seen appears markedly contrary to the impression 
engendered by. the report of your committee. It would be fair to say that manv 
wo rea. it might have thought that nuclear electrochemistry was perhaps the 

reTarTit " UniverSall y Poured to be among Jhose who do not 


'll 

^ ' 


In view of. the direct confirmation of these phenomena from so many 

b^withdrawT^ J t0 the ne § ative ERAB report on Cold Fusion should 

it be h d S an , °f flclal U ‘ S - Government document, it is important that 

it be beyond reproach (and above all, not subject to ridicule). The recent 

contenr 61 ™ 6 ^ ^ publlCations are contrary to the general drift of its 



Dr. John Huizenga 
December 26, 1989 
Page 2 


In support of this contention, relatively large amounts of money (EPRI, 

$2 M, U.S. Government Labs ~ $1 M, Japanese MITS Group c $5 M^*) anC^slated to 
be spent in 1990 in the investigation of Cold Fusion. Even the Indian 
Government has decided to fund more than 50 research workers in this area. 

The ERAB reports recommendation against large scale funding will fetter 
the U.S. competitive position in the development of a field which may, in 10- 
20 years, be seen as more important than high T superconductivity. At 
present, the Hot Fusion research budget is around $320 M (with total U.S. 
spending on fusion said to be ~ $1 B/year) . I wrote a chapter in a book on 
this subject in 1980 and recently discussed with a number of scientists active 
in this field. Both my 1980 analysis (where I showed that approach to the 
Lawson criterion had become assymptotic) and the opinions expressed to me in 
1989, suggested that no economic energy source from hot fusion is in sight. 
Those with whom I discussed the matter agreed that even if one far off day 
sustainment should be achieved, the cost of the resulting energy would be far 
more than that from fission. 

These matters, then, stress the need for a radical revision of the ERAB 
Report. The Great Problem of Fusion - sustainment - has been now achieved 
under the influence of electric fields in place of high temperatures. It is 
true that the achievement is halting and as yet unsatisfactory, but the energy 
density of the best results is that of a nuclear fission reactor - and fusion 
(i.e., tritium production) has been sustained (in occasional successful 
experiments) over many weeks. 

I put it to you that the first and necessary step is that the negative 
ERAB Report, which , in the light new research findings, no longer corresponds 
to the new findings, should be formally withdrawn. After four or five years 
of good research funding may be the time to write a new report. 

I write this letter to you not in the sense of an adversarial challenge 
but once which asks a sober and relaxed reconsideration in a subject which may 
eventually have an importance greater than that seen within DOE circles for 
hot fusion. 

With the compliments of the season, 


Sincerely, 


JOMB/jl 


J. O'M. Bockris 


Estimated from statements made to me by leading Japanese electrochemists 
in Kyoto, September 1989. It was stressed that this is only the government 
funding and that undeferred further research occurring in industry. 
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Received: by YALEVM (Mailer R2.03B) id 8531; Sun, 28 Jan 90 19:04:03 EST 
Date: Sun, 28 Jan 90 17:39:05 EST 

From: "MOSHE GAI , (203)432 5195, FAX: (203)432 3522" <GAI@YALEVM> 

Subject: Hit Probability 

To: "R.L. Garwin" <rlg2@watson> 

cc: K.G. Lynn, S.E. Jones 

J.P. Schiffer, J. Huizenga, D. Morrison. 


Dear Dick, 


Your way of thinking on hit probability is so different than mine that I took a 
hard look on it. From day one I am used to thinking in terms of the binomial 
distribution, but you avoided this by knowing that the binomial distribution at 
the limit of large N (N>7 or so) is identical to the /p 

Poisson ditribution (as well as off course to the Gaussian distribution). j 



This implies that your way of calculating hit probability should agree with 
mine only for large number of hits. Indeed for M>15 we agree within 0.8% on 
the sum of all hit probability from K=2, 3,... These are the probabilities of 
interest for the "cold fusion" experiment. \ Jr 




Out of intelctual curiosity I studied the values for small multiplicities which 
are goverened by the low fold (K) values. Indeed for multilplicity 4 you 
callculate 3.77% chance of 2 or more ring counter hits and I calculate 3.12% 

(an agreement at the level of 20%). The agreement gets better very fast, for 
multiplicity of 6, you calculate a 7.69% chance of getting 2 or more hits in 
ring counters and I calculate 7.03% (an agreement at the level of 10%). For 
M=10 K>2 the agreement is at the level of 3%. 

I find these results very fascinating, and since your way of thinking is so 
different than most of us, I thought it might be nice to write a very short 
(say a brief report) paper on the two methods, the relation between the two 
methods and a comparison of the results from the two methods. I hope you will 
find it worth your while to publish this very short and simple, but most 
fascinating paper. 


i 



l 


I will send you a first draft some time soon. Hope to hear your comments on it. 


A few other questions: 

1. May I remind you to mail to me the paper by :Menlove, Garcia and Jones for 
the NSF/EPRI meeting? (where it was stated that only 5% of the running time in 
the Los Alamos - BYU experiment was in fact relevant for the search of 
bursts). As you can imagine this fact is very important for my 
estimation of the number of bursts in our experiment, but Jones neglected to 
mention this fact to me. 


2. Did Steve Jones sent to you (as you requested) a copy of his note to 
me from last Friday where he notifies me that he will not join our publication? 
You may expect in addition a very thick package (in single space, Font 16.7) 
including a very detailed criticism of the Yale-BYU-BNL experiment. You can 
judge for yourself, however I will state that Jones and I agreed that the only 
major obstacle for publication is the value of the efficiency of our array. 

A problem now solved by your intervention. 

Some of BYU T s criticism is valid, but it amounts to correction at the level of 
a few percent. I have yet to see an experiment without problems, that can not 
be improved at the level of a few percent. Other criticism arise from not 
understanding the experimental setup. If given a chance I can demonstrate all 
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the above . 

3. I was happy to hear from Dave Lindley (Nature Washington office) that Laura 
Garwin (Nature London office) is your daughter. I had many discussions with 
her and she always left me impressed of her intelligence wit and sense of 
humor . 

Best Regards Moshe Gai. 
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Date: 25 January 1990, 17:20:39 EST 

From: (R.L.Garwin (914) 945-2555) RLG2 at YKTVMV 

IBM Fellow and Science Advisor to the Director of Research 
P.0. Box 218 
Yorktown Hts, NY 10598 
To: GAI at YALEVM 

JONESSE at BYUVAX 
KGL at APSED0FF 

Sub j ect : Conclus ions . 

Reply-To : RLG2 at WATSON 

Right now I find no significant disagreement between the efficiencies 
calculated by A1 Anderson, myself, and Moshe Gai for 2 or more votes 
among the 10 ring counters, or for 4 or more votes. 

Actually, from my tables, efficiency for 2 or more votes among the ring 
counters (with M=40 neutrons in a burst) would be 80.55% rather than the 
at least 90%" cited by Gai at the foot of page 11 (Draft 2, received by 
me 01/23/90). A1 Anderson quotes 90.45% for 50, 45.69% for 20 neutrons; 
my table shows 89.16% for 50 neutrons-- hardly a significant disagreement. 

So if people can avoid arguing about who said what when, and who said 
what first, they can get down to the significance of the results. Naturally, 
there should be a reasonable derivation of efficiency in the paper; I offer 
the view that all these numbers are small enough that there is no need for 
power-series expansions that raise a red flag of "probabilities exceeding 
unity." I would think that the paper should quote some of the Monte Carlo 
results of Anderson in comparison with analytical results. 

What does it mean? The Los Alamos results are not made more significant 
by including long-duration counting, if the bursts always come early in 
the run; so the comparison should be made with the maximum rate observed 
at LANL at the early parts of the run-- compare with burst RATE rather 
than with total bursts divided by total time. The Menlove, Garcia, Jones 
paper LAUR 89-3633 (10/16-18/89) states (Page 5) 

* ••• the majority of neutron bursts occur during warmup from LN; however 
the time period while, the sample is between -100 C and 0 C represents 
less than 5% of the sample counting time."" and 

"It is highly significant that ALL of the high-yield bursts occurred 
during the LN warm-up at -30 C." 

From Table II one can see 19 bursts with neutron numbers greater than 50- - 
many MUCH greater than 50, for which the ring counter should give 2 or 
more votes with probability exceeding 90%. Taking into account the 50% 
probability of 2 or more votes with neutron numbers as small as 22, one 
can simply say that no neutron bursts exceeding 40 occurred within the 
apparatus at Yale during the running time of the experiment. And one 
could go on to say that if the conditions were the same as at LANL, 
some xx events exceeding M = 40 would have been expected. 

Of course, if only 2 events were anticipated, even identical source 
conditions would have detected the anticipated rate with an efficiency 

l-e**-2 

since there would be a probability e**-2 that none of the two expected 
events occur. This small-number probability is much more important than 
the inefficiency of the detector. 



PRBIT SCRIPT Q1 dated 90/01/25 17:57:07 


Page 2 


Of course, if the conditions can be created only at LANL, Yale doesn’t 
have a chance... But that’s what publishing experimental results is 
all about. 

Dick Garwin 
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Received: forom CEARN by CEARN . cern . ch (Mailer R2.05) with BSMTP id 6012; Wed, 

20 Dec 89 ^4:|37:35 GVA 

Received: f/om dxmint . cern . ch by CEARN. cern. ch (IBM VM SMTP R1.2.2MX) with TCP; 

Wed, 20 DApJSJP 14:37:31 GVA 
Received: by dxmint . cern. ch (cernvax) (5.57/3.14) 

id AA11875 ; Wed, 20 Dec 89 14:22:31 +0100 
Message-Id: <8912201322 . AA11875@dxmint . cern. ch> 

Date: Wed, 20 Dec 89 14:36 GMT +1 

From : M0RRIS0N%VXPRIX . decnet . CERN@cernvax 

Subject: Cold Fusion and RF. 

To: rlg2@yktvmv 

X-Vms-To: MINT: : "rlg2@yktvmv ,T 


Dear Roberto, 20 december 1989. 

Was very interested in your account of your work and results. 

It was the first clear indication I had heard of RF perhaps explaining 
neutron bursts. As Dick Garwin had talked about such things, I sent him a 
message yesterday. He was so interested that within the hour he had phoned 
me back! He believes it is quite possible and would like to encourage you very 
much to continue as you have everything basically ready to enable you to do 
a definitive experiment. 

He also had not heard of any such experiment but told me that Menlove at 
Los Alamos had worried about such interference effects and hence had taken 
precautions to shield his experiment so that he could even put a Tesla coil 
against the outside of his apparatus and there was no effect. Thinking about 
it later, this is a negative test and shows a Tesla coil at that place has 
no effect - but it is a limited test. What you are doing is a positive test 
in that you are actually looking for RF and seeing some possible indications. 

So Menlove f s work does not necessarily rule out your conclusions - so go to it! 

Also you have a gamma counter that gives no signal whereas Menlove follows 
the Pathological Science tradition of looking for only one effect at a time 
and not investigating to see if the logical consequences (gammas, tritium, 
helium 3 and 4, are there. 

Will be in the States 1 to 21 January but will look in my Email. 

Merry Christmas and a Happy New Year 
Douglas . 
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Received: from CEARN by CEARN. cern. ch (Mailer R2.05) with BSMTP id 7128; Tue, 

19 Dec 89 18:12:54 GVA 

Received: from dxmint . cern . ch by CEARN . cern . ch (IBM VM SMTP R1.2.2MX) with TCP; 

Tue, 19 Dec 89 18:12:50 GVA 
Received: by dxmint . cern . ch (cernvax) (5.57/3.14) 

id AA00434; Tue, 19 Dec 89 17:57:54 +0100 
Mess age -Id: <8912191657 . AA00434@dxmint . cern . ch> 

Date: Tue, 19 Dec 89 18:11 GMT +1 

From : MORRI SONXVXPRIX . decnet . CERN@cernvax 

Subject: Cold Fusion. 

To: rlg2@yktvmv 

X-Vms-To: MINT: : "rlg2@yktvmv" 


Dear Dick, 19 December 1989 . 

Yesterday Roberto Battison of Perugia came to see me about their 
results . They tried to repeat the Scaramuzzi type experiments and looked for 
gammas but found that the rate of the 2.204 MeV (Bismuth) peak to the 2.224 MeV 
peak was the same whether the cell was on or off. 

Next they looked with a 3He counter and a gamma counter. They observed 
bursts of neutrons of 3 or 4 orders of magnitude above background but 
no gamma bursts at the same time. This was when Titanium shavings were 
exposed for the first time - even though the loading was very small, 

D/Ti = 0.004. Looking in more detail, they found the neutrons came in 
random spurts inside the burst. Futher when they looked at the actual charge, 
they found it was bipolar, i.e. each "pulse" consisted of a positive peak 
balanced by a negative peak so that the total area was zero. 

They were at first astonished by these results and now strongly suspect 
that there may be RF emission from the Titanium as the deuterium goes in. They 
stopped in September and have not published. 

I strongly urged them to start again and determine definitely if there is 
RF emitted. I also recommended them to adopt your suggestion of putting a 
microphone against the cathode and trying to pick up any cracks acoustically. 

It seems to me very reasonable that there could be RF emission which some 
equipment could pick up. Have you heard of anyone looking for it? If you 

could let me have any comments soon, it would be a help as if I have time 

(not at all sure), will try and write No. 21. 

Do you know if any significant changes have been made in the "draft" of the 
final report of the DOE Panel? - for example that peer-review be recommmended? 
Am a bit disappointed that the recommendations of the final report are very 
similar to those of the interim, but weaker, on the other hand it is much more 
complete and justified. 

Had hoped that Norman Ramsey would follow the tradition and after his 
Stockholm visit, pass by CERN to give a lecture, but hear his health is not 

too good - great pity. We had a lecture from Dehmalt yesterday - the first half 

was very good. 

Am going on 1st January to Houston for the DPF90 meeting which has an 
excellent programme, then to Berkeley for our Collaboration meeting (where I 
will also give my new lecture "N-Rays, Cold Fusion and Pathological Science" 
which was quite well received last week in London. Then to Utah where Steve 
Jones has invited me to BYU (he and Moshe Gai are argueing rather strongly). 
This morning received a letter from Martin Fleischmann who says he received 
all the papers that I had sent him (though does not comment on them. He says 
they have completed their big paper on calorimetry and confirmed their earlier 
conclusions. He says he will be in Lausanne in January, but I may miss him if 
I am in the States . 

Merry Christmas and All Best Wishes for 1990 (hope it is calmer than 
1989 - though this has been a good year in many ways) 

Douglas . 
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